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Thermoelastic noise will be the most significant noise source in advanced-LIGO interferometers 
with sapphire test masses. The standard plan for advanced-LIGO has optimized the optics, within 
the framework of conventional mirrors, to reduce thermoelastic noise. Recently, we and our collab- 
orators have proposed going beyond the bounds of traditional optics to increase the effective beam 
spot size and thus lower thermoelastic noise. One particular proposal for mirror shapes ( "Mexican- 
hat mirrors") yields the class of "mesa" beams. 

In this paper, we outline a general procedure for analyzing light propagating in individual arm 
cavities, and the associated thermoelastic noise, in the presence of arbitrary optics. We apply these 
procedures to study the Mexican-hat proposal. Results obtained by the techniques of this paper 
were presented elsewhere, to demonstrate that the Mexican-hat proposal for advanced-LIGO both 
significantly lowers thermoelastic noise and does not significantly complicate the function of the 
interferometer. 

PACS numbers: 04.80.Cc,07.60.Ly,42.55.-f 
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I. INTRODUCTION 

Thermoelastic noise will be the dominant noise source 
in the advanced LIGO interferometers for frequencies at 
which they are most sensitive, should the advanced LIGO 
upgrade follow the baseline design ^ Other noise 
processes (e.g., conventional bulk thermal noise; coating 
thermal noise; unified optical noise) could be important, 
but are expected to be smaller in noise power, by a factor 
~ 9 (at the minimum of the noise curve), than thermoe- 
lastic noise. As a result, the sensitivity of the advanced 
LIGO interferometers can be significantly increased. 

In a companion paper — the mesa-beam interferometer 
paper (MBI) summarized briefly in a short LIGO 
technical document (MBI-S) — we and our collabora- 
tors have argued that a signiflcant reduction in thermoe- 
lastic noise can be achieved by using modified optics that 
reshape the beam from a conventional (gaussian) profile 
into a mesa-beam shape. We furthermore performed cal- 
culations which imply (though do not comprehensively 
demonstrate) that this proposal would not place signif- 
icant burdens on the experimental community (e.g., in 
the need for slightly improved tilt control systems and 
improved mirror figure accuracies). 

The companion publications (MBI and MBI-S) survey 
the physically and practically relevant features of our 
proposal — that is, they demonstrate that the proposal 
should work, and sketch how to implement it. Those 
publications intentionally omitted many supporting de- 



tails, to be provided in thispaper and a companion pa- 
per by Erika D'Ambrosio [g]. In particular, this paper 
provides the computational details underlying our eval- 
uations of thermoelastic noise in various interferometer 
configurations; both this paper and D'Ambrosio 6] pro- 
vide computational details of our studies of the infiuence 
of mirror imperfections (tilts, displacements, figure er- 
rors) on interferometer performance. 

Further, this paper considers a broader class of possi- 
ble improvements than MBI. The MBI paper discussed 
using mesa beams reflecting off of the standard advanced 
LIGO cylindrical mirrors. In this paper, we analyze mesa 
beams reflecting off more generically shaped mirrors — 
not merely cylinders of different dimensions, but also 
frustum-shaped mirrors.^ 



A. Outline of paper 

In this paper, we provide details behind the calcula- 
tions discussed in MBI. For the reader not familiar with 
the principles behind thermoelastic noise, cavity optics, 
and cavity perturbation theory, the next four sections 
(Sections Ull IIIII Hvl and lV|) provide a brief review of the 
equations and computational techniques we employ. To 
be more explicit, in Sec^l'^G describe how to construct 
an idealized interferometer based on mesa beams while 
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A frustum is a geometric shape arising between two parallel 
planes that intersect an axisymmetric cone perpendicular to the 
cone's axis. 
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remaining within LIGO design constraints (i.e. losses per 
bounce). [In the process, we also briefly review the the- 
ory of and introduce notation for optical cavity eigen- 
fimctions and propagation operators.] We assume this 
interferometer, which has mesa beams resonating within 
the arm cavities, is driven by conventional (i.e. Gaussian) 
lasers, as described in Sec. lIIII In Sec. II VI we discuss how 
to analyze, via perturbation theory, the influence of mir- 
ror errors on interferometer performance. [In Section lVlII 
we will apply this general formalism to mirror transverse 
position errors (displacement), orientation errors (tilt), 
and figure error (shapes).] Finally, in Sec. we de- 
scribe how we numerically implemented the calculations 
described in Sections UTI and Hvl With those tools estab- 
lished, in Sections IVII and IVIII we closely follow the ar- 
guments of MBI, carefully explaining the details behind 
those arguments. 

Specifically, in Sec. IVII we use explicit forms for the 
mesa beam and its diffraction losses presented in Sections 
mand^to determine which mesa beams can operate in a 
given cavity, when the mesa beams are reflecting off of the 
standard advanced LIGO mirrrors (i.e. cylinders of ra- 
dius 15.7 cm and thickness 13 cm), as well as off mirrors 
of more generic shape (i.e. other cyliinders and frustums). 
We then employ an expression for the thermoelastic noise 
associated with a given mirror and beam configuration 
(developed generally in Sec. HTland discussed practically 
in Sec. to determine the thermoelastic noise associ- 
ated with a given configuration. Table IIVI summarizes 
our results for designs which produce lower thermoelas- 
tic noise than the baseline advanced LIGO design. 

Section IVIII describes how sensitive an interferome- 
ter using Mexican-hat mirrors will be to perturbation. 
More explicitly, Sections lyimiVII Bl and IVII Gl demon- 
strate by way of explicit comparison that arm cavities 
using spherical and mesa-beam mirrors appropriate to 
advanced LIGO are roughly equally sensitive to pertur- 
bations. Section lVlI CI describes how sensitive the overall 
interferometer (as measured by its cavity gain, its dark 
port power, and its thermoelastic noise) will be to tilt 
perturbations. Finally, Section IVII El determines how 
sensitive (as measured by cavity gain, dark port power, 
and thermoelastic noise) a mesa-beam advanced LIGO 
design without signal recycling would be to mirror figure 
errors. 

Augmented by additional results discussed in MBI, 
which demonstrate that the full interferometer including 
signal recycling will be only moderately more sensitive, 
we conclude mesa-beam interferometers could be used in 
advanced LIGO without serious concern. 



Guide to reader 

The MBI paper reports the results of calculations de- 
veloped in this paper and in one by Erika D'Ambrosio 
emphasizing their practical significance. On the other 
hand, this paper emphasizes the method by which those 



results were obtained. As a result, this paper has a dra- 
matically diff'erent structure than MBI: two of the first 
sections (Sections Ull and llV|) outline the equations that 
we have solved; the next section, Section describes 
how we solved them; and two of the last sections, Sec- 
tions!^ and [^H] briefly summarize the specific technical 
results quoted in MBI and explain how those results were 
obtained. 

The reader who wants only a cursory survey of the 
techniques used to obtain the MBI results should skim 
Sections II-IV, emphasizing the summaries in Sections 
III El and IIV El this reader should completely skip over 
the implementation section (Sec.0). Indeed, we recom- 
mend that a reader interested primarily in better un- 
derstanding results presented in MBI should work back- 
wards, flrst establishing a cross-reference between a re- 
sult in MBI and a result in Sec. IVII or IVIII of this paper, 
then reading in this paper the material surrounding the 
cross-reference. 

On the other hand, a reader who wishes to verify or 
generalize our computations should read this paper from 
beginning to end. 

B. Connection with other published work 

The MBI paper 0] and the MBLS LIGO technical 
document 0| survey our results, emphasizing their prac- 
tical significance to advanced LIGO interferometer de- 
sign. This paper provides details behind many relatively 
lengthly analytic or numerical calculations which those 
papers touched upon only briefiy. 

Another of our MBI collaborators, Erkia D'Ambrosio, 
has performed many of the same calculations (e.g., the 
sensitivity of the interferometer to tilt and to defects) 
using a sophisticated, standard, and trusted tool — the so- 
called "FFT code" — designed and developed specifically 
to study the behavior of general interferometers jQ] . As 
noted in MBI and in her paper, our two independent 
approaches agree well. 

C. Notation 

1. Symbols 

b The natural diffraction length scale (i.e. Fresnel 
length) associated with the previous two param- 
eters: b = y/L/k — 2.6 cm. 

Cv Specific heat of test mass, per unit mass, at con- 
stant volume 

D Mesa beam characteristic scale, Eq. 0. 

E Young's modulus. An elastic parameter in the 
stress tensor for an isotropic material [Eq. (|10bp ]. 
(See Appendix for a discussion of the explicit 
value used.) 
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ETM End test mass of an arm cavity (i.e. the mass op- 
posite the end Ught enters the cavity). 

Fo The net force apphed to the front face of the the 
mirror in the thermoelastic model problem of Sec. 
rrrRl cf. Eqs. © and ©. 

/ Frequency, units s^"'^. Typically, we are interested 
in values of / near the peak sensitivity of LIGO 
(i.e. / « lOOHz). 

gab Metric of 3-space. (See comments on tensor nota- 
tion, below). 

/ii_2 Height of mirror I's (or 2's, respectively) surface. 

/ Thermoelastic noise integral, Eqs. lO and 

ITM Input test mass of an arm cavity. 

k The wavevector of light in the cavity, k = 2tt/X 
with A = 1064nm. 

kb Boltzmann's constant [1.38 x lO^^^g cm^ s^^K^"^] 

L The length of the LIGO arm cavity [L — 
3.99901km]. 

£ Some number associated with diffraction losses. 
When used alone, denotes the total diffraction 
losses associated with one round trip through the 
cavity. When given certain subscripts, as with 
Ci, denotes the clipping-approximation estimate of 
diffraction losses associated with a single reflection 
off of mirror 1. 

M Mass of the test mass [40 kg] . 

MBI The mesa-beam interferometer paper 

P(r) The intensity distribution of laser light in a cavity, 
normalized to unity (i.e. J P{r)d{area) ~ 1). 

f A vector transverse to the optic axis of the cavity 
(i.e. two-dimensional). 

s Transverse displacement of an arm cavity's ETM. 

Tab Elastic stress tensor for an isotropic medium, Eq. 
(|10b|l . (See comments on tensor notation, below). 

u Electric field on some plane in the arm cavity 
(typically, the face of the ITM, propagating away 
from the ITM), renormalized to unit norm (i.e. 
/ |M|^d(area) — 1). Subscripts identify the specific 
system u refers to. 

j/° The displacement vector field associated with an 
elastic distortion. 

X, y The cartesian components of the 2- vector r. 

z A coordinate measured along the optic axis of the 
cavity (i.e. in a direction perpendicular to r). 



z is a unit vector pointing in the positive z direc- 
tion; Za = z" are the components of this vector in 
our coordinates. For example, in carteisan coordi- 
nates Za — I if a = z and if a = a;, y. 

ai Coefficient of linear thermal expansion, assumed 
isotropic. For an isotropic material, the differential 
change in volume with temperature is dV/dT/V = 
3ai. (See Appendix 0) 

ai_2 Norm of the first- and second-order corrections to 
the optical state of the cavity when the mirror is 
tilted [i.e. in Eqs. (EU and 

K Thermal conductivity, assumed isotropic (see Ap- 
pendix EJ- 

T] An eigenvalue of Eq. (|13|l . 

A An eigenvalue of Eq. . 

p Density of test mass. (See Appendix IXjl . 

a Poisson ratio. An elastic parameter in the stress 
tensor for an isotropic material [Eq. I|10b|l ] 

9 An angle through which the mirror is tilted. 

Q The expansion associated with a displacement field 
= Va2/° [e.g., in Eq. I^]. Also denotes a unit 
step function [e.g., Eq. (|llc|l ]. 

LJ = 271 f. 

Ci^2 Norm of the first- and second-order corrections to 
the optical state of the cavity when the mirror is 
displaced [i.e. in Eqs. and (|S5|l ]. 

Used only in Eq. H43f) . Cm denotes the mth zero of 
Ji{x). 

2. Tensor notation 

The elasticity calculations (e.g., involving the elastic- 
ity tensor Tab) involve tensors. While the authors pre- 
fer to interpret these expressions with an abstract tensor 
notation (i.e. Wald 0), the reader may without loss of 
generality use global cartesian coordinates and interpret 
the latin indices a and b as running from 1 ... 3, index- 
ing the coordinate directions; in these coordinates, upper 
and lower indicies are entirely equivalent (e.g. j/° = ya)- 
In this case, the reader should interpret Va as the co- 
ordinate partial derivative d/dx°' and gab as the 3x3 
identity matrix. 

3. State and operator notation 

As explained briefly in Sec. Ill C II because propaga- 
tion of the state of the field from one plane to another 
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is merely a unitary transformation, we substantially sim- 
plify the many integral-equation operations by using a 
quantum-mechanics-style notation for the operators. For 
example, the integral operation 



E{r',z')= d^r U{f',z';f,z)E{r,z) 



is represented as 



E{z') = U{z',z)E{z) 



(1) 



(2) 



Further, typically the relevant planes (i.e. the z, z') are 
known, and omitted, in expressions like the above. 

Under a similar philosophy, we also use quantum- 
mechanics-style notation for the values of the electric field 
on a specific plane; in other words, we represent states of 
the optical cavity using quantum-mechanical state nota- 
tion. For example, the optical state of the cavity may be 
equivalently denoted by the explicit functional form u or 
u{r), or alternatively as the state vector \u). In a similar 
spirit, we index the potential solutions to the eigenequa- 
tion [Eq. by an integer p, and denote the set of basis 
states as \p) . More generally we can represent any func- 
tion defined on a two-dimensional plane — not necessarily 
normalized; not necessarily a potential state — as a state 
vector. 

We use this notation in particular to simplify notation 
for inner products between two states \u) and \v): 



{v\u) — / d^r v*{r)u{r) . 



(3) 



Finally, to simplify the frequent expressions that involve 
the norm of a state, we use the symbol || • || to denote 
norms: 



(4) 



Many quantum mechanics textbooks written over the 
past half century employ this notation (cf., e.g., 0, 



II. THEORY OF AN IDEALIZED 
MEXICAN-HAT ARM CAVITY 

In this section, we review the equations needed to un- 
derstand the design of an arm cavity which (i) uses mesa 
beams, (ii) reduces thermoelastic noise in an advanced 
LIGO interferometer, and (iii) satisfies advanced LIGO 
design constraints (i.e. the diffraction losses per bounce 
are less than the advanced LIGO design threshold of 
lOppm, and the mirrors have mass 40 kg). 

Specifically, we introduce the mesa beam, defined for 
the purposes of this paper by Eq. ((SJ). We define the mir- 
ror surfaces which confine this beam by Eq. JB)). From 
Eq. (|SJ| we can determine the beam intensity profile. We 
can relate the beam intensity profile to the power spec- 
trum of thermoelastic noise [Eq. {T))] via Eqs. © and 




7.5 10 

radius, cm 

FIG. 1: The heavy solid curve is an example of a mesa beam's 
intensity profile; the lighter curve is the intensity profile of a 
Gaussian beam with similar losses. More precisely, this curve 
shows the energy fluxes per unit area (normalized to unity) 
for the mesa beam (thick curve) and Gaussian beam (thin 
curve) that have lOppm of difi^raction losses on a mirror with 
coated radius 15.7 cm. 



(|10|l . Finally, we describe how to limit attention to only 
those mesa beams with low diffraction losses, using an 
approximate approach [the clipping approximation, Eq. 
which we will test against an exact expression [Eq. 



A. Mexican-hat mirrors and mesa beams 

The mesa-beam interferometer (MBI) paper Q pre- 
sented explicit forms for certain specific mesa beams 
("canonical mesa beams") produced inside a cavity 
bounded by two identical mirrors. [In MBI, only the val- 
ues of the beam fields at the surface of the ITM (specif- 
ically, values propagating away from the ITM surface) 
were given; Appendix ^ provides a more comprehensive 
discussion.] We summarize MBI Eqs. (2.3), (2.6), and 
(2.7) to provide a compact summary of how to compute 
the values Umosa of a mesa beam's normalized electric 
field on the mirror face, at radius r from the optic axis: 



U{r,D) ^ / d^/exp 

Jr'<D 



{r-r'Y{l + i) 



2fo2 



N^{D) = / \U{D,r)\^2Tirdr 
Jo 

u„,,UD,r) = U{D,r)/N{D) 



(5a) 

(5b) 
(5c) 



where in the first integral the integration is over all points 



with r' 



< D. [For readers who wish to numerically 



explore mesa fields themselves, note that MBI Eq. (2.5) 
provides a more efficient means of calculating the mesa 
amplitude function u{D,r).] Figure ^ provides an exam- 
ple of a mesa-beam intensity distribution \u{D,r)\'^. 

As discussed in Appendix ^ and in MBI Sec. II D 
[cf. MBI Eq. (2.13)], the mirrors that refiect mesa beams 
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of the mirror; the integral is over Va, the volume of 
test mass A; and the average (...) is over one oscilla- 
tion period 1//. [Note that the dimensions of I a are 
length/force^ = s*g~^cm~-'^.] 

In the analysis that leads to the above expressions [Eqs. 
{Tj) and ©], LT 12] have made certain idealizations and 
approximations, some better than others.^ We shall em- 
ploy this expression as it stands, using the precise pa- 
rameter values provided in Appendix fXl 
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FIG. 2: The heavy solid curve is an example of a Mexican- 
hat mirror; the lighter curve is an example of a spherical (i.e. 
parabolic) mirror. More precisely, this curve shows the spher- 
ical and Mexican-hat mirrors which produce the beams de- 
scribed in Fig. 



1. Effective static elasticity model for I a 

For the frequencies at which advanced LIGO will be 
most sensitive (/ ~ lOOHz), the interior of the test mass 
will respond effectively instantaneously to the imposed 
pressure profile in the model problem described above. 
As a result, by going to the accelerating frame of the test 
mass and following that accelerating frame for one period 
to simplify the average, we can approximate I a by 



back into themselves (denoted Mexican-hat mirrors) nec- 
essarily have a continuous height function hyui given by 



(6) 



where k — 27r/A is the wavenumber of light in the arm 
cavity. The solid curve in Figure |21 provides an example 
of a Mexican-hat mirror. 



B. Thermoelastic noise 

Using the expressions in the previous section, one can 
design idealized interferometers which operate with mesa 
beams and Mexican-hat mirrors. The thermoelastic noise 
associated with the resulting interferometer design can be 
discovered with the techniques of Liu and Thorne (LT) 
|T^. who employ Levin's direct method |0| to derive the 
following formula for an interferometers' thermoelastic 
noise in terms of a noise integral I a'- 



Shif) = AnkB 



aiET 



2 4 



Cv(l - 2a)pujL 



(7) 



cf. LT Eqs. (3), (4) and (13); the notation used here has 
been described in Sec. II CI 

The sum on A = 1, 2, 3, 4 is over the four test masses, 
and the quantity summed is the noise integral 



I A = 



{VQf^ dvolume 



(8) 



In Eq. ©, 6 is the expansion (i.e. fractional volume 
change) resulting from a pressure — sinusoidally oscillat- 
ing at frequency / — a;/27r, proportional to the light 
beam intensity, with net force Fq — imposed on the face 



Ia 



(Ve)^dvolume 



(9) 



Va 



where the expansion 8 — Va?/° arises due to the re- 
sponse (i.e. local displacement field) y° of the test-mass 
substrate to the following static conditions: 

1. Static (accelerating-frame) force: Some net force 
FoZ acts uniformly throughout the test mass (i.e. 
provides a force density Fo/Va)- To be specific, 
using the equations for an isotropic elastic solid (cf. 
Blandford and Thorne 14]), the response field ?/° 
satisfies 



Tab = 



E 



(l + cr)(l -2(7) 

(1 - 2a) 



(10a) 



(10b) 



For example, they ideahze the test-mass material as an isotropic 
material; however, the proposed test mass material, sapphire, 
is only moderately well approximated as isotropic. They also 
assume the elastic response of the solid is quasistatic (introducing 
errors of order the sound crossing time over the gravitational 
wave period, ~ lO""^) and that heat diffuses slowly so the elastic 
response of the solid to the imposed pressure profile is adiabatic 
(introducing errors of order the diffusion length over the beam 
size, or ~ 0.3 mm/10 cm ~ 3 X 10^'^). 

The mirrors considered in this paper have characteristic dimen- 
sion of order H ^ 10 cm. Thus, the sound crossing time is of or- 
der H/cs ~ H/^E/p ~ 10-5 sec. Since the body responds elas- 
tically to imposed forces on times of order a few sound-crossing 
times, the elastic response to a force which is imposed at fre- 
quency ~ lOOHz is effectively instantaneous. Equivalently, the 
same conclusion follows because the gravitational wave frequency 
is far below the frequency of any resonance of the test mass. 
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2. Static pressure on mirror face: An equal and oppo- 
site net force —F^z acts on the mirror face, with a 
distribution P{r) proportional to the beam inten- 
sity profile. As a result, the displacement field y°- 
must satisfy 



Tab{r,z^O)z' - -FoP{r)za , 



(10c) 



where the mirror's front (reflecting) surface is at 
z = and where the mirror lies below z = 0. 

3. Break translation symmetry: Finally, translation 
symmetry must be broken (i.e. so numerical sim- 
ulations converge to a unique solution). To break 
translation symmetry, we fix "the location" of the 
mirror, where that location is determined as the 
average location of points in a set TZ. Therefore, to 
specify a unique solution, we require the displace- 
ment field satisfy 



(lOd) 



where TZ is some arbitrary nonempty set and 
denotes the average over TZ. 



In 



[This quasistatic approximation to I a and the overall 
thermoelastic power spectrum [Eq. iQ] involves errors 
only of order 10~^ relative to the general Liu-Thorne 
expression Q for 1^ 0|; cf. footnotes [3 and O or Ap- 
pendix |0] 



Propagation, the eigencondition, and diffraction 
losses 



In the previous two sections, we described mesa beams 
and the technique by which one can determine, given a 
beam and four mirrors that support that beam (two iden- 
tical ITMs and two identical ETMs), the thermoelastic 
noise an interferometer using those mirrors and beams 
would experience. But not all combinations of mesa 
beams and mirrors satisfy advanced LIGO design con- 
straints; the advanced LIGO design specifications require 
each arm cavity have low diffraction losses per bounce 
(lOppm). 

In this section, we describe how the diffraction losses 
can be precisely computed [Eq. (|15|l ] and estimated [the 
clipping approximation, Eq. (|16|) ]. Using the diffraction 
losses, one can then find combinations of mesa beams and 
mirrors which produce low diffraction losses per bounce. 

Our analysis proceeds by describing in general terms 
the propagation of light in an arm cavity, focusing in 
particular on eigenmodes and their loss per round trip. In 
the process of justifying our simple estimate of diffraction 
loss, the clipping approximation [Eq. we will review 

the definitions and properties of eigenmodes of the arm 
cavity. This review establishes notation and conventions 
for our presentation of optical perturbation theory, as 
discussed in Sec. IIVI 



1. Principles of paraxial optics 

The light in LIGO arm cavities is well-described by 
the paraxial approximation. In this section, primarily 
to establish notation conventions, we briefiy review the 
propagation of light under the paraxial approximation. A 
detailed description of the relevant physics can be found 
in standard references.^ Briefly, in this approximation, 
we can completely describe the state of the wave by the 
wavelength A of the light used and the values of the wave 
amplitude on some fixed plane z ^constant (i.e. trans- 
verse to the optic axis) and at some time t (cf. 14] ) . The 
values at any causally-related later combination t' , z' — 
which must be separated from the transverse plane z 
at time t by a light ray (i.e. t' — t = \z — z'\/c, mod 
reflections) — can be deduced by applying the appropri- 
ate linear functional to these states. 

For example, we could characterize the optical state 
by the value of the electric field on some plane z.^ We 
would then relate the field at any other plane and at any 
other time to our initial state via the linear operation 

E{r',z',t')^ I d^rU{r=',z';f,z)E{f,z,t), 



cf. Eq. iQ. For compactness and clarity, we employ a 
quantuni-mechanics-motivated notation,^ in which the 
integrals are suppressed and the operation above is de- 
noted by [cf. Eq. 101] 

E{z',t') = U{z',z)E{z,t) . 

By way of example, the following are kernels of integral 
operators which describe free propagation down a length 
L] reflection off a mirror of height ft. 1,2; and a "window" 
that cuts out all light outside a radius r = -Ri,2, respec- 
tively: 



GL(f,r') 

Gi,2(r,r') 
ri,2(r,r') 



■ exp I 



(f. 



+ kL 



'2ttL ' [ 2L/k 
— 5(r — r ') exp [— 2ifcft,i^2(?')] 
S{r^r')e{Ri.2-r) 



(11a) 

(lib) 
(11c) 



where Q{x) is a step function equal to 1 when x > and 
otherwise, and where the negative sign in Gi arises 
because of boundary conditions on the electric fleld at 
a perfectly reflecting surface; cf. 0,0,0, and other 



^ For a more pedagogical and yet brie f pre sentation, we recommend 
Blandford and Thome's treatment 1141 . Other pedagogical (and 
tech nica l ) treatmen ts can be found in many laser physics books, 
e.g., ll5lll6llT^ll8l . We also recommend the collection of origi- 
nal rese arch articles on laser and resonator physics compiled by 
Barnes Il9l . 

^ We limit attention to a single polarization, e.g., the polarization 
excited by LIGO. 

^ Most quantum mechanics textbooks written over the past half- 
century (e.g. Slol llCl ) have adopted a similar operator notation. 
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GLG2GLk> G2GLk> 




FIG. 3: This figure demonstrates how propagation operators 
described in Eq. IIH are used to evolve the field from one 
plane of constant z to another. For brevity, the operators 
T\ and T2, which allow for the finite extent of the mirrors' 
surfaces, have been omitted in this figure. 



the operator Gh\ These eigenmodes are resonant when 
j^g^ikL j,g^^ g^j^^ positive. The arm cavity length L and 
wavevector k are tuned so one state, the ground state Uo, 
is resonant; in this case, rjo, k, and L satisfy 



Arg [Tyoe^''^^] = . 



(14) 



The mesa beam [Eq. (0)] is designed specifically to be 
one of the eigenmodes of the Mexican-hat mirrors [Eq. 
©] when the mirrors are infinite (i.e. Ti = T2 = 1), 
and when the cavity length is properly tuned to admit 
it. When the mirrors are finite, the true eigenmodes of 
an arm cavity (be it bounded by spherical or Mexican- 
hat mirrors) will have a slightly different shape than the 
shape supported by infinite mirrors. Further, the eigen- 
value of the arm cavity 77 will no longer be of unit mag- 
nitude (i.e. \r]\ < 1). 



previously noted references on paraxial optics for further 
details. 

These propagation operations can be combined to gen- 
erate more complicated processes. For example, while 
Gi describes reflection off an infinite mirror of height hi, 
GiTi describes the reflection of light off a finite mirror 
of height hi and radius Ri. 

2. Describing paraxial propagation through an arm cavity 

When we have an arm cavity bounded by two cylin- 
drical mirrors separated by a length L, of cylinder radius 
Ri^2 and surface heights hi^2 respectively, we can use the 
above propagators [Eq. Hll|l ] to describe the free, un- 
driven propagation of light within the arm cavity''. For 
example, as we demonstrate graphically in Fig. |31 we can 
express the field at the surface of mirror 1 — specifically, 
the field for light heading away from the mirror surface — 
at time t + 2L/c in terms of the field at that surface at 
time t as follows: 

u{t + 2L/c)^GiTiGLG2T2GLu{t) (12) 

for the operators on the right side as given in Eq. Hll|l and 
for \u) denoting some polarization of the electric field. 

We are in particular interested in eigenmodes — that is, 
states u so the beam returns proportional to itself: 

T]e'^'''^u = G1T1GLG2T2GLU (13) 

for some rj. [We factor the phase e^*'^^ out of ry to 
eliminate the effect of the phase factor e*'^^ present in 



The description of an arm cavity driven by a source laser adds 
only a straightforward inhomogeneous term to our discussion; 
the eigenmodes of the homogeneous term are as always in such 
problems of paramount importance. 



3. Diffraction losses: Exactly 

From Eqs. H12|l and H13|) . we know that when the ini- 
tial optical state is an eigenmode, the optical state after 
n round trips decays as measured by its norm, which 
evolves as \\u{t + n2L/c)\\ = |7y|"||u(i)||. Thus, ff |r;| < 1 
the magnitude of u (i.e. the norm of u) decreases with 
each round trip. 

Since our model above only permits losses from diffrac- 
tion, and since the power in the cavity is proportional to 
/ d2r|u|2 (i.e. the 

norm of u, squared), the following 
quantity is the diffraction loss per round trip: 

Cnet = 1 - . (15) 

When the two mirrors are identical, as is our case, we 
can subdivide the loss in two, to obtain a meaningful 
diffraction loss "per bounce" . 

We can therefore extract from numerical solutions to 
Eq. p3|l (cf. Sec. IV B(l the magnitude of 77 and, by 
the above expression [Eq. IjlSI) ]. deduce the mesa-beam 
diffraction losses in the presence of finite mirrors. We do 
not, however, use this method in this paper. 

4- Diffraction losses: The clipping approximation 

In practice the exact expression for the diffraction 
losses given above [Eq. l(TB|) ] is slow to evaluate, because 
we must solve for fine details of the eigenequation for 
each mirror size i?i_2 of interest. To obtain a rough esti- 
mate of the diffraction losses, usually accurate to a factor 
of order unity, we will often use instead the clipping ap- 
proximation. 

The clipping approximation estimates diffraction losses 
by assuming the beam is not significantly changed by 
those losses; the beam has to a good approximation the 
same shape as it would have in the presence of infinite 
mirrors. Assuming the beam profile is known, we then 
directly compute the losses by determining the power lost 
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FIG. 4: A contour plot of the (logj^Q of the) diffraction losses 
of the mesa beam [Eqs. IBlll with z = -L/2, i.e. IIB3ll 1 accord- 
ing to the clipping approximation [Eq. I|16|l ]. as a function of 
mirror radius R and averaging scale D. Solid curves indicate 
losses 10"" for n = 2, 3,... 11. The two heavy solid curves 
correspond to n = 2 and n = 5; these two curves have special 
significance to LIGO design (they determine the maximum 
value of losses off the back and front faces of an ITM, per 
bounce; cf. Sec. Ill Dl . 

off the edges of finite mirrors on each reflection. Specifi- 
cally, we assert £nct ~ -Cciip, with 

-Cciip = >Ci + >C2 (16a) 
£i = |(l-Ti)«„|2 = I \uJ^Sr (16b) 

Jr>Bi_ 

C2 = \{1-T2)GlUo\^ = f \GLUo\^d\ il6c) 

Jr>Ri 

where Uq is an eigensolution to the propagation equation 
[Eq. (|13|l ] computed assuming infinite mirrors. 

5. Diffraction losses of the ideal mesa beam, computed via 
the clipping approximation and corrected to estimate the true 
diffraction losses 

In this paper, we approximate the diffraction losses 
of the mesa beam resonating between two identical 
Mexican-hat mirrors by the clipping approximation [Eq. 
(|16|l ] applied to the mesa beam [Eq. ©]. These losses 
are shown in Fig. 

The true diffraction losses per bounce generally dif- 
fer from the estimates of clipping approximation. For 
example, for & D = Ah mesa beam resonating between 
-R = 16 cm cylindrical mirrors, the true diffraction losses 
are about 18.5 ppm, while the clipping approximation 



estimates diffraction loses of 21 ppm. Based on this 
good agreement, in this paper we assumed the true 
diffraction losses were well-approximated by the clipping- 
approximation estimate, with relative errors of order a 
factor 1.2. 



D. Advanced LIGO design constraints 

The advanced LIGO interferometers are composites of 
many interrelated systems; each component of that sys- 
tem has been designed making assumptions about the 
other components. In particular, the advanced LIGO in- 
terferometers require (roughly) the following constraints 
on mirror and beam designs: 

1. Mirror mass o/40kg: We consider only 40kg test- 
mass mirrors, because the present advanced LIGO 
design for the suspension and seismic isolation sys- 
tem requires 40kg mirrors. [Also, the advanced 
LIGO suspension system cannot support heavier 
mirrors (Phil Willems, private communication)].* 

2. Mirror front radius limited by arm- cavity diffrac- 
tion losses: We consider only combinations of mir- 
ror front face radius R and beam size (e.g. if mesa 
beams are used, D) which have diffraction losses 
equal^ to 10 parts per million (ppm).^° 

3. Mirror back radius limited by power-recycling-cavity 
diffraction losses: Finally, we consider test-mass 
mirrors whose back face is sufhciently large that 
light entering the arm cavity through the back face 
loses at most 1 percent (i.e. 10, OOOppm) of the in- 
put power. ^"'^ When evaluating this constraint, we 
assume the input light is in the same state as the 
arm cavity light (i.e. in a cavity eigenmode). 



The support system could hold up lighter mirrors. However, 
thermoelastic noise decreases as the mirrors and beams are in- 
creased in size (in equal proportions). We therefore limit atten- 
tion to the most massive possible mirrors. 

The advanced LIGO design could accept mirrors with lower 
losses. However, since thermoelastic noise decreases (as a gen- 
eral rule) when the beam radius increases, we limited attention 
to the largest possible beams. 

This constraint is required to keep the diffraction losses bounded 
by a reasonable portion of the total loss. Specifically, for the 
baseline design there is 125 W of input power to the interferom- 
eter and 830 kW of circulating power in each arm cavity. The 10 
ppm of diffraction loss per bounce results in a diffraction power 
loss in the arm cavities of 4 X lOppm X 830kW = 33 W, which is 
25 per cent of the 125 W of input light, a reasonable value. 
As with the constraint involving the front face size (cf. foot- 
note ^21: this constraint is required to keep the diffraction losses 
bounded by a reasonable proportion of the total loss. Specifi- 
cally, for the baseline design the power impinging on each ITM 
is 1.05 kW, so one per cent diffraction loss (i.e. 10, OOOppm) cor- 
responds to losing 2 X 0.01 X 1.05kW = 21 W at the ITM input, 
which is 17 per cent of the 125 W total laser power. 
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[This constraint only matters for noncylindrical 
mirrors (i.e. frustums) with a relatively small back 
face size; for cylindrical mirrors, it holds automat- 
ically.] 



E. Summary: Exploring mesa-beam arm cavity 
configurations 

To summarize, we design and evaluate new mesa-beam 
interferometer configurations by the following process: 

1. Pick a configuration: We (i) select two axisym- 
metric mirrors (i.e. cylinders or frustums) which 
have the same front face radius and (ii) select some 
value D. We shape the mirror faces [Eq. @] and 
tune the length of the arm cavity [Eq. H14|l ] so the 
ground state eigenmode of that cavity — the eigen- 
mode which is approximately a mesa beam with 
length scale D — is resonant. 



III. DRIVING A MESA-BEAM ARM CAVITY 
WITH A GAUSSIAN BEAM 

Ideally, a mesa-beam arm cavity would be driven by 
light already in the mesa-beam state. But if the optics 
required to generate mesa-beam inputs are unavailable 
or too inconvenient, Gaussian beams — generated by con- 
ventional optics — can be used to drive the arm cavities. 
With the proper choice of Gaussian, the arm cavity be- 
haves almost as if the ideal mesa beam was used. 

In the remainder of this paper, as in the mesa-beam 
interferometer paper (MBI) jj| , we shall assume properly- 
chosen Gaussian beams are used to drive the interferom- 
eter. 

In her companion paper |^, Erika D'Ambrosio pro- 
vides a comprehensive analytic and numerical discussion 
of Gaussian beams driving mesa-beam cavities [cf., e.g., 
ED'A Eq. (3)]. In this section, we briefly summarize some 
of her results on mesa-beam cavities driven by gaussian 
beams. 



2. Check that the configuration satisfies the advanced 
LIGO design constraints: If (i) the mirrors masses 
are not equal to 40 kg, if (ii) the physical diffraction 
losses in the arm cavity are greater than 10 ppm per 
bounce [Eq. ifT^ ]. or if (iii) the diffraction losses 
for input light are greater than one percent, then 
we stop and try again: the configuration does not 
satisfy the advanced LIGO design constraints de- 
scribed in Sec. Ill Dl [In this paper, we use only 
the clipping approximation, Eq. (|16|l . to estimate 
diffraction losses.] 

3. Evaluate the thermoelastic noise of this configura- 
tion: We finally evaluate the noise integrals Ii and 
I2 for each of the two arm-cavity mirrors [Eqs. © 
and (|10|l ]. using as beam pressure profile the (nor- 
malized) beam intensity profile [Eq. lO]:^^ 

P{r)^W,,s.[r,D)\^ . (17) 

[Section fV Al will describe in detail precisely how we 
evaluated the two thermoelastic noise integrals.] 

These thermoelastic noise integrals then determine 
[via Eq. Q] the overall thermoelastic noise of an in- 
terferometer using two identical arm cavities, where 
each cavity consists of this particular mirror and 
beam configuration. 



A. Beam and cavity used in this section 

1. Fiducial mesa-beam cavity 

Rather than study all possible interferometer configu- 
rations (i.e. all possible mirror and beam sizes), we select 
a single fiducial mesa-beam arm cavity: the cylindrical 
mirrors have front face size 16 cm and the mesa beams 
have D = 4b = 10.4 cm, where b = ^XLjl-K = 2.60 cm 
(and where A = 1.064/im and L = 4km are the light 
wavelength and arm length) [cf. MBI Sec. IV A 2]. This 
fiducial arm cavity has diffraction loss 18 ppm.^^ 

We also assume the mirrors bounding the arm cavity — 
both the input test mass (ITM) and end test mass 
(ETM) — have no intrinsic losses, except for diffraction 
losses (which we treat separately). Specifically, we as- 
sume the ETM is perfectly reflecting, and assume the 
ITM has power transmissivity t\ = 0.005 [ij. Thus, our 
power reflection and transmission coefficients for the ITM 
and ETM are 

r\ = 0.995 , t\ = 0.005 ; 

r| = 1.0 , i| = 0. (18) 

[This same fiducial cavity will appear frequently else- 
where; for example, we explore perturbations of this par- 
ticular fiducial cavity when exploring the sensitivity of 
the interferometer to tilt.] 



^ Ideally, we should use as beam intensity profile the true eigen- 
state appropriate to two finite mirrors, as obtained by solving Eq. 
1131 . However, because diffraction losses are low, the eigenstate 
for finite mirrors is very well approximated by the eigenstate 
for infinite mirrors; and the latter is the mesa-beam amplitude 

^mcsa ■ 



While not quite appropriate — the diffraction losses are slightly 
larger than the advanced LIGO design cutoff of 10 ppm y| — 
this fiducial design is sufficiently similar to acceptable advanced 
LIGO designs. 
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2. Gaussian beam used to drive the fiducial arm cavity 

As described in MBI Sec. IV B, this fiducial arm cavity 
will be driven by a Gaussian beam of form given in MBI 
Eq. (2.8), with Gaussian beam radius = 6.92cm. As 
in MBI, we denote this state by \ud)- This state has 
overlap 



(19) 



with the arm cavity's ground state Uo [cf. MBI Eq. (4.2)]. 
[Erika D'Ambrosio has found that this particular Gaus- 
sian beam gives the largest possible overlap with our fidu- 
cial beam; cf. ED 'A Eq. (3).] 

Strictly speaking, we use a Gaussian state \ud) which 
is shifted by precisely the phase offset necessary to make 
7o real and positive; therefore. 



7o 



VO.940 = 0.970 . 



(20) 



B. 



Operation of fiducial cavity when driven by a 
Gaussian beam 



Equation (|19|l implies a fraction 0.94 of the (gaussian) 
driving-beam light will enter the Mexican-hat-mirrored 
cavity; the remaining 6 percent will be reflected. The 
fraction of the light that enters the cavity is amplified 
(due to resonant interaction with the bounding cavity 
mirrors), but eventually leaks back out the same way 
it entered (minus that fraction lost through diffraction 
losses). The combination of the 6 percent reflected light 
and 94 percent transmitted light (modulo losses) then 
returns back towards the beamsplitter. 

This small section provides a more quantitative view 
of this entire process. 



1. Decomposing light along the cavity eigenstates 

The transmission and reflection processes are best un- 
derstood in terms of the resonant eigenstate of the cavity. 
We therefore rewrite the driving beam in terms of a pro- 
jection along and perpendicular to the ground state: 



\Ud) = Vo\ud) + il-Vo)\ud) 

^ Jo\Uo) + il-Vo)\ud) 



(21) 



where the projection operator Vo is defined so VoV 
\uo) (uolv) for any state v, i.e. 



Vo = \Uo) {Uo\ , 

and where 70 is defined by (cf. Sec. IIII A 2(1 

7o = {uo\ud) ; 



(22) 



(23) 



the phase of Ud has been designed so the unperturbed 70 
is real. Since the input state Ud is normalized. 



\VoUd\\^ + \\il-Vo)ud\\^ = \\Ud\\' 



1 



(24) 



The portion joUo of the driving field excites the mesa- 
beam cavity, while the portion (1 — Vo)ud cannot res- 
onate and thus gets fully reflected off the input test mass 
(ITM). 



2. Amplification of light entering cavity 

A fraction 7^ of the light power enters the arm cavity, 
resonates, and is amplified by a factor 



l + ri 

l-n 



1 



2Co 

l-n 



(25) 



within the arm cavity, where Co denotes the diffraction 
losses per bounce for the ground state of the arm cavity. 
Therefore, relative to the overall input power, the light 
in the arm cavity is amplified by a factor 



^anal To 



1 -frj 



1 - 



2Lo 



737. 



(26) 



[While the term £o/(l — r/) is very small in this case and 
can be ignored, we provide the entire expression now for 
use in Sections IIVI and IVIII when we study the effect of 
mirror perturbations on the cavity gain.] 



3. Output: Transmission of light out of cavity, and 
recombination with reflected light 

The cavity's output light is the combination of light 
refiected at the ITM (which does not enter the cavity) 
and light exiting the arm cavity through the ITM. After 
some algebra (cf. 0),^^ we find that given a driving beam 
Ud decomposed according to Eq. H21|) the refiected light 
Ur is in the state 



\Ur) = Vo\ud)-{l-Vo)\ud) 

= {2Vo - 1) \ud) ■ 



(27) 



where for simplicity we have ignored the (small) effects of 
diffraction losses on the light that resonates in the cavity. 



The second (unnumbered) equation in D'Ambrosio gives the 
transmitted field through the arm cavity. Dividing this by the 
second mirror transmissivity, squaring it to get an expression for 
beam power, and then speciaUzing to r2 = 1 and r\ = rj, we 
obtain the first factor in Eq. I25i . The second factor arises when 
losses are properly taken into account on each bounce. 
The portion of the beam amplitude proportional to the arm cav- 
ity ground state enters the arm cavity and resonates. It therefore 
exits that the arm cavity with the same phase it enters. On the 
other hand, the rest of the light experiences nearly perfect de- 
structive interference. For the rest of the light, the arm cavity 
acts like a perfectly reflecting mirror, so that portion of the light 
amplitude exits the arm cavity with the opposite phase than it 
entered with. For a more precise demonstration of this effect, 
see for example the first equation in D'Ambrosio with r2 = 1, 
t\ = tj and ri = rj, and both exponential factors unity. 
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IV. EFFECT OF MIRROR PERTURBATIONS 
ON ARM CAVITIES AND THE 
INTERFEROMETER 

In Section we described how to design an idealized 
mesa-beam arm cavity and to evaluate the thermoelastic 
noise associated with such a design. 

In this section, we shall describe in extremely general 
terms how one can use perturbation theory to study the 
influence of defects on a single ideal (i.e. lossless) arm 
cavity. We make no reference to the specific details of 
the cavity used, save that — for technical convenience — 
we limit attention in our perturbation theory to a system 
with two identical, infinite mirrors. [Recall the canon- 
ical mesa beams, as defined in Sec. Ill Al are designed 
for symmetric cavities.] Naturally the methods for per- 
turbing arbitrary mirrors bear considerable similarity to 
methods for perturbing purely spherical mirrors. Many 
explicit tools and concepts carry over with little change 
to our more generic approach. 

Specifically, in Sec. lIV Al we introduce the orthonormal 
basis set of solutions we shall use to construct the per- 
turbation theory expansion. Then, in Sec. IIV Bl we list 
the explicit second-order perturbation expansions we use 
to deduce the effect of mirror defects upon the optical 
state of the cavity. Finally, in Sec. IIV CI we describe how 
changes in the optical state of the cavity influence other 
quantities that can be deduced from that state, such as 
the cavity diffraction losses [Eq. (|16(l ] and the thermoe- 
lastic noise integrals I a [Eq. • 



A. Basis states for perturbation theory 

We will use perturbation theory to analyze the effect 
of small changes in the arm cavity mirrors on the solu- 
tions to the eigenequation ((T^ . To construct perturba- 
tion theory expansions by methods similar to those used 
in quantum mechanics (cf., e.g., iSt lEe iO] ) > we prefer to 
use as a basis for the perturbation expansion the eigen- 
modes of the initial equation [Eq. (|13|l ].^^ Fortunately, 
as shown in Appendix^ so long as the two mirrors are 
identical and inflnite, the set of eigenmodes of Eq. l(T^ 
are a complete, orthogonal set, independent of the mirror 
shape. 

At this point, we could introduce basis states (and no- 
tation for basis states) which are defined as eigensolutions 
to Eq. H13() . For technical reasons — our methods simplify 
practical evaluation of the perturbation equation^^ — we 



prefer instead to use a basis designed to simplify "half" 
(i.e. a square root of) the round-trip operator, with a 
certain phase factor removed. In other words, we use as 
basis states the eigensolutions \p), Xp to the equation 

Xp \p) = e-^^-^GiGL \p) (28) 

where \p) is some state and p denotes some index over 
all basis states. [By an argument following that given in 
Appendix^ this equation admits a complete set of or- 
thogonal solutions \p).] These states correspond directly 
to solutions to the full eigenequation. Explicitly, if we 
insert a solution \p) to Eq. H28|) into the eigenequation 
[Eq. H13|l ]. we immediately conclude that \p) is also an 
eigensolution of Eq. (|13|l . with eigenvalue 

Vp = Xl. (29) 



B. Effect of perturbations on light in the arm 
cavities 

When the mirror shapes hi^2 are deformed, the light 
propagating in the cavity changes. Given a basis of states 
and a specific problem to perturb [Eq. I|13|l ]. we employ 
conventional techniques from quantum mechanics (cf., 
e.g., 0,0,1^1) to compute that change. 



1. Results of perturbation expansion 

Using conventional quantum-mechanics-style tech- 
niques (cf., e.g., lESEj) but note the operator we per- 
turb is unitary rather than hermitian), we can devise a 
perturbative expansion for the cavity ground state eigen- 
value and state u'^ in powers of the height perturbation 
62- The derivations of the sometimes-long expressions 
noted here are provided in Appendix lEl 

For simplicity, we provide perturbation theory expan- 
sions only to second order, and only to describe the ef- 
fects due to changes in the height of one mirror (i.e. /12 of 
mirror 2) by an amount dh. To express these changes in 
height in dimensionless form, we introduce the variable 
52: 

62 = 2k6h2 . (30) 
Change in eigenvalue: The eigenvalue r]o of the ground 



^® Actually, as Erika D'Ambrosio has frequently reminded us, we 
can construct perturbation theory just as well in the more general 
case, by using a dual basis. See, for example. Appendix A of 
D'Ambrosio I3. 

^'^ On the one hand, by removing the common factor exp(ifcL) 
present in Gl, we produce a naturally dimensionless eigenvalue 
problem (as k and L enter only through the single length pa- 



rameter b). On the other hand, by focusing on "half" of the 
propagator rather than a round-trip operator, we can easily de- 
duce the relationship between the fields at the two ends of the 
cavity. The latter proves helpful, because we will study perturba- 
tions of the mirror at one end of the cavity (end 2) , yet represent 
the state of the cavity field by its values on the other end (end 
!)■ 
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state \uo) = |0) changes as 



fe#0 



Vo - Vk 



(31) 

Change in eigenstate: When we construct the pertur- 
bation theory expansion to second order, we find the 
ground state changes according to an expression of the 
form 



o,pt/ 



)«|0) + |v^W) + |^/.(2)^ . (32a) 



In this expression, tp^^^ and i/'^^^ denote those terms first- 
and second-order in 52 in the perturbation expansion, re- 
spectively. To be exphcit, when we perform the pertur- 
bation theory expansion [details of which are provided in 
Appendix IE] . we find 



(1) 



fe#0 



Vo - Vk 



\k){k\52\0) , (32b) 



(2) 



o'^k 



Vo - Vk 



^ 'lo 



Vk 



(fc|<52|0) (0 1^210)-^ 



"Tip 



p/0 



Vo - Vp 



{k\62\p){p\S2\0) 



(32c) 



r 



By construction jV''-^-') and jV'''^'') are orthogonal to the 
unperturbed ground state |0). As a result, this expression 
(|32|1 for If/)') is not normalized: we find, working to second 
order, that the norm of p^) is 



«ptl<pt)«l + llV'^'^lP + (5('52') 



(33) 



where we use the shorthand = Therefore, 

the physically appropriate normalized perturbed state 
is given by the expression 



1 



(1) 



(2; 



|o) 



(1) 



(2) 



|0) 



(34) 
0{5l). 



2. Estimating convergence of perturbation expansion 

Perturbation theory is only effective when higher or- 
der terms provide only a small correction to lower order 
terms. To test the convergence of the series, we compare 
the magnitudes of the first two perturbative corrections. 
When < 1, 

we believe the series converges 
and our expressions should be effective. 



C. Implications of change in optical state for other 
quantities 

In the previous section, we described the effect of per- 
turbations 5hi^2 on light propagating in a single Fabry- 
Perot arm. These perturbations cause the beam shape 
incident on the two mirrors to change, generally in a dif- 
ferent way at each mirror. Therefore, quantities that 
depend on the state of the beam at each mirror, such as 
the diffraction losses and thermoelastic noise integrals /a. 



also change. In this section, we loosely describe how the 
lowest order effect of these changes can be characterized. 



1. Change in diffraction losses under perturbation of one 
mirror 

Unfortunately, a systematic treatment of diffraction 
losses within the context of perturbation theory proves 
very tricky, not the least because we must represent both 
the height change and the effect of finite mirror size as 
perturbations, then use many states to insure the pertur- 
bation theory expansion converges properly for the effects 
of diffraction. 

Rather than perform a truly accurate, well-motivated 
computation, we will in this paper loosely estimate 
diffraction losses by the clipping approximation applied 
to the perturbed beam state. By way of example, we can 
estimate the diffraction losses produced during a single 
reflection off of mirror 1 or off mirror 2 in the presence of 
a perturbation of mirror 2 by expanding the appropriate 
clipping approximation estimate Ci or C2 [cf. Eq. Hl()|) ]: 



£1 « {u\Oi\u) , 

£2 « {Glu\02\Glu) 

01 = 1-Ti, 

02 ^ I-T2. 



(35a) 
(35b) 
(35c) 
(35d) 



To be explicit, we can expand the clipping-approximation 
estimate for £1 to second order in 82 as 

A « (0|Oi|0)+2Rc((v(i)|Oi|0)) (36) 

+ (^W|Oi|^w)-||v«ino|Oi|o) 

+ 2Re((v^'^|Oi|0)) +0{5l) . 
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A similar expression is found for £2, with Oi replaced 
by G\02Gl- [Because the beam profile will generally 
change in different ways at the two ends due to a pertur- 
bation localized at only one end, generally Ci 7^ £2-] 



2. Change in cavity gain due to perturbations of one mirror 

The amount of power resonating in the cavity also 
changes, in part through modified diffraction losses [Eq. 
(|35|l ] but also through a change in the resonant ground 
state of the cavity from Uq (i-e. |0)) to u'^ [Eq. Be- 
cause the resonant state of the cavity changes, the overlap 
7o changes to 7^, which to second order can be approxi- 
mated by 



7o = {u'o\ud) = 
where we define 



1)||2 



7o+7i||^('^ll+72||^^'^|| , 
(37a) 



7i = (i/^(i)|u<i)/||V^'^ll 



(37b) 
(37c) 



To be very specific, we can find the perturbed cavity 
gain S^^jjjji by evaluating Eq. to second order in Sh2, 



l-ri 



1 - 



i-n 



(38) 



using Eq. for the perturbed diffraction losses £1 and 
£2 and using Eq. ifTfl) for the change in overlap between 
the perturbed cavity ground state and the driving beam. 

[Though we do not write out the resulting second-order 
expansion in full detail, we will apply these methods to 
compute to second order the effects of mirror tilt on cav- 
ity gain in Sec. IVffI ] 



3. Change in dark port power due to perturbations of one 



An interferometer consists of two arms. When the 
beamsplitter recombines the two fields leaving the two 
arm cavities, the light going out the dark port of the in- 
terferometer (udp) is the interference between the light 
reflected off the two arms (denoted I and 11) : 



(39) 



If both cavities are identical, then m^ji — Urj = Ur for Uj. 
given by Eq. (|28|l and no light exits the dark port. If the 
cavities are perturbed, however, power will generically go 
out the dark port 



For example, if only cavity 11 is perturbed, then Urji 
u' , = Ur, and the dark port beam state is 



1 



"dp 



{Uj, - Ur) 



V2' 

= V2(7X-7oMo) 

= V2 [(-7„| 1^^(^)112+ 71 II^W II +72||V'(2) \uo) 



(70 +71 II V^^'^ II 



(1) 



+70 



0(^2') 



(40a) 



where going from the first line to the second we use Eq. 

[for Ur]; from the second to the third we use Eqs. 

[for Po] and ^ [for 70]; and from the third to the 
fourth we use the perturbation expansions (|32|) [for u'^] 
and ^ [for 7^,]. 

The dark port power Pdp (as a relative fraction of in- 
terferometer input power, which is twice the input power 
to each arm cavity) can be expressed as 

^dp = ll"dp||V2, (40b) 
an expression which we shall not expand here. 



Influence of perturbations of one mirror on 
thermoelastic noise 



The thermoelastic noise integrals associated with each 
mirror (1 and 2) in a given arm cavity both change be- 
cause the beam profile at each mirror changes. To evalu- 
ate those changes, we insert the modified state [Eq. H34|l ] 
into the thermoelastic noise integral Ia [Eq. 0]. 

For example, to evaluate the linear-order change in 
thermoelastic noise at mirror 1, we first take the new 
normalized intensity profile P{r) at mirror 1, given by 



Piir) ^ \u'o{r)f 



\u{r)\'+2Re(u*{r)i;^'Hr))+0{dl) .(41) 



Using this new intensity profile, we solve the thermoe- 
lastic noise model problem [Eq. (|10|) ] for the expansion 
O' and in particular the first-order change in expansion 
SQ = Q' — Q. Finally, we insert this first order change SQ 
into the expression for 1^ with A = 1, linearized about 
the background intensity profile: 

6Ii = ^ [ (Vae)(V"(5e)civolume (42) 
Fq Jvi 

where O is the expansion produced by the unperturbed 
pressure profile FqP. 

The beam shape at mirror 2 also changes; by a similar 
construction, we can find its effect on the thermoelastic 
noise integral l2- 
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D. Observations which simplify our computations 

1. Symmetry and the influence of first- order changes in 

state 

The ground-state beam of physical interest is isotropic. 
As a result, when we perform the computations outlined 
above, we find only the axisymmetric (about the optic 
axis) part of I^A*-^-*) contributes to first-order changes in 
the three integral quantities of physical interest (the over- 
lap 7o; diffraction losses; and thermoelastic noise). 

By way of example, consider the first-order change in 
thermoelastic noise [Eq. H42|l ]. Since the unperturbed 
pressure profile and thus unperturbed expansion (0) are 
isotropic, the above integral couples only to the axisym- 
metric part of 59 and therefore SP and ■0'^^ [cf. Eq. H41|l ]. 

Similarly, to lowest order displacement and tilt per- 
turbations have odd parity; therefore, symmetry insures 
that the first-order change 71 [Eq. (|37|l ] is zero. 

2. Assuming (roughly) equal contributions from the two 

mirrors to changes in I a and diffraction losses 

The change in mirror shape at mirror 2 produces 
roughly comparable changes to the beam profile at both 
ends. To a rough approximation, then, we can assume 
that the changes in thermoelastic noise and diffraction 
loss at mirror 2 will be comparable to their changes at 
mirror 1. 



E. Summary: Using perturbation theory to 
explore the sensitivity of the arm cavity and 
interferometer to perturbations 

To summarize, we evaluate how sensitive mesa-beam 
arm cavities and interferometers are to mirror perturba- 
tions as follows: 

1. Assume the arm cavity is fiducial: The perturba- 
tion expansions above can be applied to any arbi- 
trary mirror configuration (i.e. to any specific mesa- 
beam D)}^ In principle, we could apply pertur- 
bation theory to all possible mirror configurations 
generated in Sec. ^1 Instead, to avoid repeating 
computations which should yield nearly identical 
results, we apply perturbation theory only to the 
fiducial advanced LIGO arm cavity presented in 
Section nil A II 



Our technique assumes the cavity eigenmodes of an arm cav- 
ity with finite mirrors are well-approximated by the modes of 
a cavity with infinite mirrors. Therefore, the only relevant pa- 
rameter remaining is the mesa-beam averaging scale D. [This 
approximation ignores quantities of order the diffraction losses, 
10 ppm ~ 10-5.] 



2. Consider only perturbations of one ETM: Similarly, 
we change only one end test mass in one arm cavity, 
rather than apply perturbations to each test-mass 
mirror. 

3. Compute the eigenmodes of the fiducial arm cav- 
ity: We find the natural eigenmodes of the fiducial 
arm cavity — and therefore the basis states \p) in 
our perturbation expansion — by numerically solv- 
ing the basis-state eigenequation H28|l . [The numer- 
ical code which solved this integral eigenequation is 
described in Sec. lVBI ] 

4. Find how the ground state of arm cavity changes: 
Next, we apply the perturbation expansion [Eq. 
(j^ J to find how the ground state of the arm cavity 
changes when each perturbation of physical interest 
is applied: tilt, displacement, and mirror figure er- 
ror. In other words, for each of these perturbations 
Sh2, we use Eq. (|32|l to find the first- and second- 
order corrections j^/'^"'^'') and to the state of 
the cavity. 

5. Determine how the thermoelastic noise integral Ia 
for each mirror changes: Given the changed beam 
state, we can recompute the thermoelastic noise 
integral Ia for each mirror (cf. Sections ^ and 
IV A|) . using the perturbed beam state u'g{r) and the 
perturbed beam intensity profile P'{r) = \u'^{r)\'^. 
More directly, the series expansion of thermoelas- 
tic noise can be discovered by a series expansion 
of the thermoelastic noise integral, as sketched (to 
first order) in Sec. IIV C 41 

6. Describe how the cavity gain and dark port power 
change: Finally, given the changed beam state, we 
can also recompute the arm cavity gain [Eq. H38|l ] 
and interferometer dark port power [Eq. H4UI) ] when 
the ETM of one cavity is perturbed. These two 
expressions both depend on how the perturbed arm 
cavity interacts with the driving beam [i.e. on 71 
and 72; cf. Eq. H37|) ] and on the perturbed arm 
cavity state itself [i.e. on | [V'''"'"-' 1 1 and [[■0(^^[[; cf. 
Eq. (|32|l ]. Also, the cavity gain depends on how 
the diffraction losses of the arm cavity change [Eq. 



V. NUMERICAL IMPLEMENTATIONS OF THE 
THERMOELASTIC NOISE INTEGRAL, THE 
OPTICAL EIGENEQUATION, AND OPTICAL 
PERTURBATION THEORY 

In the previous two sections, we described the abstract 
expressions we must evaluate to design and analyze an 
arm cavity bounded by Mexican-hat mirrors. In this sec- 
tion, we describe how we implemented and solved those 
equations. 
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Where independent methods were used to perform a 
particular computation, we indicate the different tech- 
niques used. 

A. Thermoelastic noise for perfect (undeformed) 
mirrors 

To evaluate the thermoelastic noise power spectrum 
[Eq. |(7I)], we need to perform two tasks. First, we must 
solve the elasticity problem described in Eq. I|10(l : then, 
using the result, we can evaluate the integral which 
enters directly in Eq. Q. 

We employed three methods to address these two tasks. 
The first two methods — a numerical (finite-element) so- 
lution and an exact analytic solution — applied only to 
special circumstances: axisymmetric beam profiles, and 
cylindrical mirrors. The third method was an approxi- 
mation based on assuming the mirror to be half-infinite; 
we used it only as a quick, easy-to-evaluate check on the 
qualitative behavior of the previous two methods. 



To be explicit, LT Eq. (44) gives an expression that 
is precisely 1/2 of / (1/2 because of averaging which we 
have factored out but which LT retain) , in terms of quan- 
tities defined in LT Eqs. (35) and (36). To generalize to 
a generic axisymmetric pressure profile, one need only 
change LT Eq. (36) — the only equation which involves 
the specific pressure profile — so that it involves the in- 
tensity function P{r) defined in Eq. ()10c|l for m > 0:^'^ 

Pni = 2 7-2//- / P{r)Jo{CTnr/a)rdr . (43) 

Here a is the radius of the cylindrical mirror (denoted 
R elsewhere in this paper); is the mth zero of Ji{x); 
and the functions Jq and Ji are the zeroth and first order 
cylindrical Bessel functions. The sum converges rapidly: 
typically, only a handful of terms in the infinite sum [LT 
Eq. (45)] are required. 

3. Approximate technique: half-infinite mirrors 



1. General approach 

We usually solved the elasticity equations HlObp - ljlOdl) 
using a commercial two-dimensional finite-element code 
p^ . We chose the region TZ in Eq. (jlOd|l to be one of 
the mirror faces. [To use a two dimensional code, we 
limited attention to axisymmetric mirrors and pressure 
profiles.] 

The commercial finite-element code we employed gave 
us the displacement vector We then used postpro- 
cessing code discussed in Appendix [H] to evaluate the 
derivatives and integrals needed in Eq. @. 



2. Special case: cylindrical mirror 

Independently, we employed the analytic elasticity so- 
lution Liu and Thorne (LT) ^3] developed for cylindrical 
mirrors with axisymmetric pressure profiles imposed on 
them. 

Liu and Thorne constructed a solution to the elasticity 
equations for a cylinder [Eqs. IjlObp - ljlOdjl ] which they ap- 
ply to Eq. to find an explicit expression for the ther- 
moelastic noise produced by gaussian beams on cylin- 
drical mirrors. By replacing a single equation in their 
expressions, we can convert their solution to one appro- 
priate to arbitrary axisymmetric beam intensity profiles. 



The code works faster if the region 7^ is a point or set of points. 
However, the code is significantly more susceptible (on physical 
grounds — the points act like "nails" in the mirror) to small errors 
in the neighborhood of points, errors that contribute significant 
erroneous expansion. Therefore, as a practical compromise we 
chose 7^ to be a surface. 



If the mirror is sufficiently large compared to the im- 
posed pressure profile P{r), the elasticity problem we 
must solve [Eq. HlU|l ] can be well approximated by a so- 
lution to a similar problem with the mirror boundary 
taken to infinity. In this C8jSG, cLS the force density term 
[Eq. H10c|l ] goes to zero, we need only solve for the re- 
sponse of a half-infinite (i.e. filling the region z < 0) elas- 
tic medium to an imposed surface pressure. As described 
in greater detail in Appendix^ analytic expressions ex- 
ist for the response, permitting us to find an compact 
expression for the thermoelastic noise integral. We find 



(l + g)(l-2g) 
2nE 



K 



P I 



(^)f (44) 



where P{K) is the two-dimensional Fourier transform of 
P{r): 



P{K) EE j d'^fe^'^-^P{f) . (45) 
Recall that P{r) is normalized to unity (cf. Sec. Ill B|l . 



B. Numerically solving for the resonant optical 
eigenmodes of a cavity bounded by arbitrary 
axisymmetric mirrors 

To test the validity of the clipping approximation [Eq. 
(|16|l ] against the exact diffraction losses [Eq. (|15|) ] and 
to generate the set of basis solutions needed to construct 



The LT coefficient po is independent of the pressure profile shape; 
it is always 
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perturbation theory [Eq. H28() ]. we must numerically solve 
for the cavity eigenmodes. In this section, we describe 
numerically how we converted the cavity eigenmode in- 
tegral equation into a numerical eigenproblem, which we 
then solved with standard numerical tools (e.g., Mathe- 
matica) . 

1. Setting up the problem to he solved; preliminary analytic 
simplifications 

Rather than solve the eigenequation for a full round 
trip through a symmetric cavity [Eq. (|13|l with i?i — 
-R2], we exploit symmetry and instead study the closely 
related problem of half a round trip through a symmetric 
cavity. In other words, we plan to diagonalize the basis 
state eigenequation [Eq. when that eigenequation 

is restricted to the space of functions defined on Ri = 
R2 = b Xmax: 

A \u) = TiGiGlTi exp{-ikL) \u) (46) 

If we write out the appropriate eigenequation for this 
"half-a-round-trip" operator, making use of the defini- 
tions in Eq. Hll|l and using the dimcnsionless spatial units 
X — r/b, we find the following integral equation: 

AV'(X) = J d^X'g{X;X')i;{X') , (47a) 

giX;X') ^ -^eiX„,,^-\X'\)e[ ' ^ Vb) 
h = 2kh. (47c) 

Here, tp{X) — u{r)/b denotes a dimcnsionless represen- 
tation of the state u. 

Since the mirror surfaces are (ideally) axisymmetric, 
the operators above all commute with rotation around to 
the optic axis. We can therefore mutually diagonalize this 
operator and the generator of rotations. Therefore, we 
require ■0 be proportional to exp(zTO(/9) for some integer 
to: 

V'(X, (^) EE $(X)e™^ . (48) 

Substituting in this form for ip^ we reduce the problem 
to a series of one-dimensional integral equations, one for 
each |to|: 




XdXg„,{X,X)<^{X) , (49a) 



= -z™+ij„,(XX)e4^^-''W](49b) 

2. Method of numerical solution 

For each to of interest, we represent the integral op- 
erator on the right side of Eq. (|49|) as a matrix. For 



simplicity, we discretize space in a uniform grid Xa = 
AX^^^/{N - 1) for A = 0... TV, define $a = H^a) 
and similarly, and approximate the integral by a simple 
quadrature rule.^^ Thus, we approximate the integral 
equation of Eq. (|49() by the matrix eigenvalue problem 

This equation can be solved for 77 and $ by any standard 
eigensolution package. 

3. Interpreting and applying the numerical solution 

To summarize, to find numerical approximations to 
eigensolutions '0 of Eq. H46|l — an eigenequation similar to 
the basis-state eigenequation [Eq. but for mirrors of 
finite radius R — bX-a^a.^ — we construct and solve the 
matrix eigenvalue problem (|50|l . 

The approximate eigensolutions rjp and $ so obtained 
provide numerical approximations to the two cavity- 
eigenproblems of technical interest: 

• Solutions for true cavity eigenmodes: This method 
provides us with precisely the numerical cavity 
modes needed to understand mesa beams in the 
presence of finite mirrors. [Specifically, eigenmodes 
and eigenvalues of Eq. are also eigenmodes of 
Eq. Ijl^ll .] In particular, we can use the norm \ri\ to 
determine the round-trip diffraction losses accord- 
ing to Eq. (113 . 

• Solutions for the basis states: Further, by making 
the mirrors sufficiently large that diffraction effects 
can be ignored'^^, solutions to Eq. 1)46(1 provide good 
approximations to solutions to the basis eigenvalue 
problem [Eq. fl^]. 

C. Numerical implementation of perturbation 
theory 

To use perturbation theory exactly [i.e. Eqs. (|31|l and 
(|32|l ]. one needs an infinite collection of states. In prac- 
tice, we limited attention to a handful: the lowest-lying 
ten to twenty for each |to| = 0, ... 7; more modes were 
used when the rate of convergence of the perturbation 
theory expansion suggested more were needed. 



We in fact used the equation Eq. l5Ui as stated. More sophisti- 
cated quadrature techniques, such as Gaussian quadrature, offer 
greater accuracy with fewer points (and hence significantly less 
computation time). 

As a practical matter, we test the quality of our solutions by 
observing the convergence of the norm \r]\ as we increase the 
number of points used to represent the state. 
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D. Numerical exploration of changing diffraction 
losses 

Given the (5/i2-induced perturbed beam state at 
mirror 1, we can compute the (cUpping approximation 
estimate of the) diffraction losses at mirror 1 for a given 
height perturbation Sh2 at mirror 2. Largely, we simply 
evaluated the integrals (i.e. inner products and norms) 
required to construct Eq. However, to provide an 

independent numerical check (i.e. to insure we had no ty- 
pographical or structural errors), we also evaluated the 
diffraction losses directly, using the definition Eq. p6|l . 
for a sequence of height perturbations 5h{e) = e5h2 with 
e = 0,0.1,0.2, . . .1; then fitted a second-degree polyno- 
mial to the resulting data points to extract the series 
coefficients in Eq. 



E. Numerical investigation of changes in 
thermoelastic noise when one mirror shape is 
perturbed 

Finally, we can use the known form of ip'^^'^ (r, ip) [from 
Eq. (|32b|) ] in the procedure outlined earlier [cf. Eqs. H41|l 
and H42I) ] to evaluate the first-order change in thermoe- 
lastic noise associated with mirror 1 due to changes in 
the shape of mirror 2. 

Specifically, given -(/^'^'(r, 0), we use Eq. 141|) to find 
how the pressure profile changes. We then use the first- 
order change 5P to the pressure profile in the LT expres- 
sion for thermoelastic noise (see Sec. IV A"2)l . linearized 
about the response to the unperturbed cavity beam in- 
tensity P. [We will not provide the rather lengthy but 
straightforward linearization of the LT expressions here; 
any computer-algebra system can easily reproduce the 
desired expansion.] 



VI. DEPENDENCE OF THERMOELASTIC 
NOISE ON MIRROR AND BEAM SHAPE 

In this section, we explore the dependence of the ther- 
moelastic noise integral I a [Eq. ©] on the various arm 
cavity parameters available to us: (i) the mirror's di- 
mensions (i.e. cylinder height and radius) and shape (i.e. 
frustum or cylinder), (ii) the beam size (e.g., the mesa- 
beam size D), and (iii) the type of beam resonating in 
the arm cavity (i.e. Gaussian or mesa). The thermoelas- 
tic noise integral provides a simple way to characterize 
how the thermoelastic noise power spectrum Sh [Eq. iQ] 
of these configurations compares to Sh for the baseline 
advanced LIGO configuration {Sf^): 

Sh/Sf^ = ///bl , (51) 

where /bl (evaluated below) is the value of the ther- 
moelastic noise integral for the baseline advanced LIGO 
configuration [Eq. H52|) ]. and where / is the value of the 



thermoelastic noise integral for all four (identical) mir- 
rors in the interferometer. 

Specifically, in this section we compare the following 
types of mirror and beam configurations, all of which sat- 
isfy the advanced LIGO design constraints [Scction lll Dj : 
(i) the baseline advanced LIGO configuration (which uses 
Gaussian beams and cylindrical mirrors; (ii) an improved 
baseline configuration (which also uses Gaussian beams 
and cylindrical mirrors); (iii) configurations with mesa- 
beam light resonating between identical cylindrical mir- 
rors; and (iv) configurations with mesa-beam light res- 
onating between identical frustum mirrors. [In this sec- 
tion, we restrict attention to arm cavities with identical 
mirrors; in Appendix^ we discuss generalizations to arm 
cavities bounded by nonidentical mirrors.] In Table IIVI 
we summarize the optimal (i.e. lowest value of I/Ibl) 
configurations we found for each class. The results of 
items tabulated in this section are applied, in MBI Sec. 
Ill (cf. MBI Table I), to produce advanced LIGO designs 
with lower thermoelastic noise than the baseline design. 

Our evaluations have been performed independently by 
all three co-authors (RO'S, SS, and SV), using multiple 
methods (both finite-element solutions and infinite-sum 
analytic solutions) when appropriate. 



A. Baseline advanced LIGO configuration 

The baseline design of an advanced LIGO interfer- 
ometer [l| has four identical cylindrical sapphire test 
masses (i.e. physical radius Rp — 15.7 cm; thickness 
H = 13 cm; mass 40 kg) whose surfaces are coated over 
most of their surface (i.e. out to a radius R = Rp—S mm). 
These mirrors' surfaces are designed so the largest possi- 
ble Gaussian consistent with the 10 ppm diffraction loss 
constraint (Sec. Ill D|l resonates in the arm cavity [i.e. a 
Gaussian beam with radius Tq = 4.23 cm — 1.636; cf. Ap- 
pendixIFland MBI Eq. (2.8)] for a discussion of Gaussian 
beams] . 

The thermoelastic noise integral for an arm cavity 
bounded by cylindrical mirrors and using a Gaussian 
beam is found by (i) constructing the Gaussian ampli- 
tude function UG{r, To) [MBI Eq. (2.8)] and its associated 
beam intensity profile P{r) = |ug(t', t'o)|^; (ii) solving for 
the elastic expansion O that arises due to P{r) in the elas- 
tic model problem of Eq. H10|l described in Sec. IIIBI (cf. 
Sec. IV AI for numerical methods); and (iii) inserting the 
resulting expansion into the definition of the thermoelas- 
tic noise integral I [Eq. Q]. For the baseline beam and 
test mass, the resulting value of the noise integral / is 

/bl ^ 2.57 X 10-^Vg-^cm-i . (52) 

[The advanced LIGO cylindrical mirrors and Gaussian 
beams are optimal: these beams produce very nearly the 
lowest thermoelastic noise possible using Gaussian beams 
refiecting off identical 40 kg cylindrical mirrors, where 
those cylinders are coated out to a radius R = Rp — 8mm 
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TABLE L The thermoelastic integral I for a cylindrical test 
mass and a Gaussian beam, in units of the value Jbl = 2.57 x 
10~^^s''g~^cm~^ for the advanced LIGO baseline design. The 
values of //7bl are estimated to be accurate to within one per 
cent. 



R 


-Rp[cm] 


H[cm\ 


To 


ro/b J/Jbl -Co 


[ppm]» 


BL** Rp 


- 8mm 15.7 


13. 


4.23 


1.63 1.000 


10 


Rp 


15.7 


13. 


4.49 


1.73 0.856 


10 



"Diffraction losses in each bounce off the test mass, in ppm (parts 
per million), computed in the clipping approximation 
'Baseline design for Advanced LIGO interferometers 



and where the diffraction losses per bounce are restricted 
to less than 10 ppm.] 



B. Improved baseline advanced LIGO configuration 

The conventional baseline described above wastes the 
last Smrn of mirror face size. We can improve upon the 
thermoelastic noise merely by eliminating the uncoated 
ring in the last 8mm, i.e. by coating the mirror out to 
the edge and adjusting the beamspot size to fill in the 
extra space. If the mirror coating extends out to the test- 
mass edge so i? = Rp = 15.7 cm, and the Gaussian beam 
radius is correspondingly increased to Vo — 4.49 cm so the 
diffraction losses are still 10 ppm, then the thermoelastic 
noise is reduced to I/Ibl = 0.856; see Tabled 

[Again, the same advanced LIGO cylindrical mirrors 
(i.e. with unchanged phyisical radius and thickness) pro- 
duce very nearly the lowest possible thermoelastic noise, 
among all 40 kg cylindrical mirrors coated out to their 
physical radius.] 



TABLE II: The thermoelastic integral / for a cylindrical test 
mass and a mesa beam, in units of the value 7bl = 2.57 x 
10~^^s''g~^cm~^ for the Advanced LIGO baseline design. The 
values of 7/7bl are estimated to be accurate to within one per 
cent. The first four test masses, like the baseline, are mirror 
coated only out to R = Rp — 0.8cm; the last five are coated 
all the way out to the test-mass edge, R — Rp. 



i?[cm 


1 


Rp [cm] 


i7[cm] 


D/b 


///bl 


-Co [PP 


Rp - 


8mm 


14.67 


14.79 


3.00 


0.414 


10 


Rp- 


8mm 


15.70 


13.00 


3.43 


0.364 


10 


Rp- 


8mm 


17.11 


10.87 


4.00 


0.442 


10 


Rp- 


8mm 


19.58 


8.30 


5.00 


1.000 


10 


Rp 




13.94 


16.38 


3.00 


0.373 


10 


Rp 




15.70 


13.00 


3.73 


0.290 


10 


Rp 




16.37 


11.88 


4.00 


0.313 


10 


Rp 




18.85 


8.96 


5.00 


0.628 


10 


Rp 




21.36 


6.98 


6.00 


1.69 


10 



"Diffraction losses in each bounce off the test mass, in ppm (parts 
per million), computed in the clipping approximation 

ther out to (i) their physical radius Rp or (ii) only out 
to Rp — 8 mm. Table Ull summarizes our results for each 
one-parameter family. 

Whether the mirror is coated out to the full physi- 
cal mirror radius or not, in both cases quadratic fits to 
I/Ibl{Rp) give minima rather near the baseline phys- 
ical radii Rp = 15.7 cm. Indeed, to within our ac- 
curacy of computation, the same mirror shape used as 
the baseline advanced LIGO design {Rp = 15.70cm and 
H = 13.00cm) gives the optimal thermoelastic noise for 
mesa beams with 10 ppm diffraction loss. The beam radii 
D and thermoelastic noise I/Ibl for these near-optimal 
mexican-hat test masses are shown in Table HVl below. 



C. Mesa beams reflecting off^ identical 40kg 
cylindrical mirrors 

Arm cavities with mesa beams reflecting off cylindrical 
mirrors admit configurations with even lower thermoe- 
lastic noise than the improved baseline. To explore the 
advantages of mesa beams, we evaluated the thermoelas- 
tic noise integral (via the method described in Sec. IIIE() 
for two one-parameter^^ families of cylindrical mirrors 
and mesa beams: 40 kg mirrors with the largest possi- 
ble mesa beams resonating off their front faces (i.e. set 
by lOppm diffraction losses; cf. Sec. Ill G'5)| .^'^ coated ei- 



The arm cavity has three free parameters (mirror radius; mirror 
thickness; and mesa-beam radius parameter D) and two con- 
straints (mirror mass and diffraction losses per bounce). 
As expected, thermoelastic noise decreased with increasing mesa 
beam radius; we obtain the lowest value of thermoelastic noise 
integral / when the mesa beam radius is as large as possible, 
consistent with the diffraction constraint. 



D. Mesa beams reflecting ofl' identical 40kg 
frustum mirrors 

Roughly speaking, two mirrors with a larger front face 
radius permit a wider beam to resonate in the arm cavity 
and yield even lower thermoelastic noise. Therefore, we 
explore arm cavities bounded by frustum-shaped mirrors 
(cL footnote^l, which expand the front face of the mirror 
at the expense of the back face. 

More specifically, we considered arm cavities resonat- 
ing with mesa beams of beam radius parameter D 
bounded by 40 kg test-mass mirrors whose front and back 
faces were as small as diffraction losses would permit (i.e. 
the front face produces precisely 10 ppm diffraction losses 
with arm-cavity mesa beam light of scale D; the back 
face produces precisely one percent diffraction losses for 
the same mesa-beam input light; cL Sec. IIID]) .^^ These 



While in principle we could consider any combination of front 
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TABLE III: The thermoelastic integral / for a frustum input 
test mass (ITM) and a Mesa beam, in units of Jbl = 2.57 x 
10~^*s*g~^cm~^. The values of ///bl are estimated to be 
accurate to within one per cent. 
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Rp2\ 


cm] 


H[cm] 
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I/Ibl 




£o[ppm] 
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10 & 10** 
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12 
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13 


.22 


13.45 


4.13 


0.208 


10 & 10'' 


Rp 


— 8mm 


19, 


.58 


15 


.27 


10.43 


5.00 


0.285 


10 & 10'' 


Rp 




13, 


.94 


9 


.88 


22.24 


3.00 


0.345 


10 & 10" 


Rp 




15, 


.70 


11 


.56 


17.00 


3.73 


0.198 


10 & 10" 


Rp 




16 


.37 


12 


.19 


15.49 


4.00 


0.175 


10 & 10" 


Rp 




17, 


.29 


13 


.04 


13.75 


4.39 


0.162 


10 & 10" 


Rp 




18, 


.85 


14 


.58 


11.33 


5.00 


0.193 


10 & 10" 


Rp 




21, 


.36 


17 


.00 


8.62 


6.00 


0.398 


10 & 10" 



"Radius R of the mirror-coated portion of both the inner and the 
outer faces of the test mass, in units of its physical radius Rp. 

'Diffraction losses in each bounce off the test mass, in ppm (parts 
per million); the first number is for the light inside the arm cavity, 
on face 1 of the test mass (radius Ri); the second number is for 
the light impinging from the beam splitter onto face 2 of the test 
mass (radius J?2)- 



TABLE IV: Optimized test-mass and light beam configura- 
tions, their thermoelastic noise compared to the baseline. [A 
subset of this table appears as Table I in MBI i.] 



Test Masses 

{fipi, Rp2; H} 




Beam Shape 


\ J TE 


BL: cylinders, R — Rp - 


- 8mm 


BL: Gaussian 




{15.7, 15.7; 13.0} 




ra ^ 4.23cm 


1.000 


BL: cylinders, R. — Rp - 


- 8mm 


mesa 




{15.7, 15.7; 13.0} 




D/b ^ 3.73 


0.364 


identical frustums, R. — 


Rp — 8mm 


mesa 




{17.11, 12.88, 14.06} 




D/b ^ 4.00 


0.207 


BL: cylinders, R. — Rp 




Gaussian 




{15.7, 15.7; 13.0} 




To — 4.49cm 


0.856 


BL: cylinders, R. — Rp 




mesa 




{15.7, 15.7; 13.0} 




D/b ^ 3.73 


0.290 


identical frustums, R — 


Rp 


mesa 




{17.29, 13.04, 13.75} 




D/b ^ 4.39 


0.162 



arm cavities satisfy advanced LIGO design constraints 
(cf. Sec. lTTDll . 

Table Hnl sumniarizes the thermoelastic integrals ///bl 
for two one-parameter family of designs: (i) the mirror 
is coated out to its physical radius R = Rp and (ii) the 
mirror is only coated out to i? = Rp — 8mm. By fitting a 



and back face radii Ri and i?2 and any mesa beam radius D 
such that all three satisfy the LIGO design constraints presented 
in Sec. Ill Dl we found that for R2 greater than or equal to the 
minimum radius allowed by diffraction losses the thermoelastic 
noise integral increases with R2 (for fixed mirror mass and front 
face size). Therefore, we limited attention to R2 as small as 
possible. Also, as usual, we limited attention to mesa beams as 
large as diffraction losses on the mirror front face permit. 



quadratic to (/)/6, ///bl), we estimate the optimal mir- 
ror dimensions and associated beam radii D. Our opti- 
mal results appear in Table Hvl 



VII. INTERFEROMETER SENSITIVITY TO 
MIRROR PERTURBATIONS 

In Sec. IVII we found mirror and beam configurations 
for the advanced LIGO arm cavity which indeed pos- 
sess lower thermoelastic noise than the baseline advanced 
LIGO design. But our primary modification — the change 
to Mexican-hat optics for the cavity arms — involves em- 
ploying mirrors which have never before been used in an 
interferometer. Naturally, then, we must make every ef- 
fort to demonstrate that this radical proposal will not 
introduce new problems. 

For example, the Mexican-hat mirror has a very flat 
central region (cf. Fig. In our early presentations of 
this proposal, it was suggested to us that such a mirror 
design might make the resulting interferometer substan- 
tially more susceptible to errors, be they from static tilts 
and displacements or mirror figure error. 

In this section, we examine this concern by examining 
the effect of perturbations on both gaussian and mesa- 
beam arm cavities. More explicitly, we by applying the 
tools described in Sec. lIVI to two fiducial beams (cf. MBI 
Sec. IV A 1): 

• Fiducial mesa beam and fiducial Mexican-hat mir- 
rors: The fiducial mesa-beam arm cavity, described 
in Sec. IIII A II which has mesa beams with D = 
46 = 10.4 cm. 

• Fiducial Gaussian beam and fiducial spherical mir- 
rors: A fiducial Gaussian beam arm cavity, which 
has beam radius Tq — 4.70 cm (i.e. a 5-value g — 
0.952).^^ [This fiducial Gaussian-beam arm cav- 
ity differs from the basehne (cf. MBI Sec. IV A 1).] 
The fiducial Gaussian beam cavity is chosen to have 
the same diffraction losses, on a mirror of the same 
coated radius, as the mesa-beam arm cavity (i.e. so 
the two fiducial cavities we compare are similar). 

Using these two fiducial beams, we demonstrate that 
mesa-beam and Gaussian-beam interferometer designs 
for advanced LIGO will have broadly (i.e. within a factor 
a few) similar sensitivity to perturbations. In short, 
we demonstrate that the mesa-beam proposal will not 
introduce undue sensitivity of the arm cavities and inter- 
ferometers to mirror errors. 



The value g = 1 — L/7?.j , for TZj the radius of curvature of 
spherical mirrors, can be related to the Gaussian beam radius ro 
using formulae presented in Appendix IFl 
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A. Frequency distribution of parasitic modes as a 
measure of arm cavity sensitivity to perturbations 

Given the denominators present in Eqs. and (|32|l . 
a system will generically be more unstable to perturba- 
tions if the eigenphases of excited modes [i.e. a,rg{r]i^)] 
are close to the eigenphase of the resonant state. We 
can therefore crudely characterize the influence of generic 
perturbations by the distribution of eigenmodes nearby 
the ground state, also called the parasitic mode distribu- 
tion. This discussion provides the basis for MBI Sec. IV 
B. 

The resonant frequencies u = kc are determined when, 
after one round trip through the arm cavity, light in a 
given state interferes constructively with itself. There- 
fore, the light must be in an eigenstate \p) of the cavity 
[Eq. H13|l ] and, moreover, the frequency of the light must 
be chosen so the eigenvalue of the round-trip eigenequa- 
tion is real, or chosen so [cf. Eq. (|14|l ] 



TABLE V: For a LIGO arm cavity with fiducial Mexican hat 
mirrors {D /b = 4, R — 16 cm), this table gives the separation 
Auj/uiFSR of the eigenfrequencies of parasitic modes from the 
eigenfrequency of the fundamental mesa-beam mode. [The 
table is indexed by two numbers, p and I, because the parasitic 
modes may be written in the form u{r) = Rp^i{r) exp{±iW); 
cf. Sec. IV Bl The index p is the number of radial nodes in the 
eigenfunction Rp{r).] 







I = 


I = 1 


I = 2 


I = 3 


p = 





0.0 


0.0404 


0.1068 


0.1943 


p = 


1 


0.1614 


0.2816 


0.4077 


-0.4581 


p = 


2 


0.4303 


-0.4140 


-0.2570" 


-0.0812" 


p = 


3 


-0.2330" 


-0.0488" 


0.1406" 


a 



27rn — 2Luj/c + Arg(r/p) 



(53) 



"If the cavity length is adjusted so this mode is resonant, then its 
diffraction losses C in each bounce off the test mass exceed 10000 
ppm (i.e. one percent, in contrast with the fundamental mode's 18 
ppm), so it cannot resonate strongly. 



B. Effect of displacement on cavities bounded by 
spherical and Mexican-hat mirrors 



for n some integer. The same eigenmode \p) resonates 
at a specific frequency, and every other frequency sep- 
arated from that frequency by the free spectral range 
t^FSR = TTc/L. Within each free spectral range, different 
eigenmodes (i.e. \p), \q)) are occur at different frequen- 
cies, separated by an amount uniquely determined by 
their eigenvalues 77: 



Aujpg ^ujp-ujq = — [Arg(?7p) 



Arg(77,)] . (54) 



For Gaussian beams, the eigenmodes are distributed 
within each free spectral range regularly; each nearest 
neighbor is separated by frequency [cf. Eq. (|1''4() ] 



Aw 



L 



9 = WFSR X 



9 



(55) 



0.099 X WFSR 



using the g-value g = 0.952 for our Gaussian baseline 
beams. 

For the mesa beams, the frequency distribution of par- 
asitic modes must be obtained numerically, by (i) solv- 
ing the eigenequation for the eigenvalues rjp [Eq. (|13|l . 
or equivalently Eq. (|28|l . using the numerical methods 
of Sec. IV B| and then by (ii) using those eigenvalues in 
Eq. if^H) to deduce Awpo and therefore the distribution of 
parasitic modes. Table fVl lists the values for Awpo/^FSR 
for a few states. Among states with low diffraction losses 
(i.e. with £ = 1 — |77p|^ < 10~^), most modes are very 
well separated from the ground state; the nearest para- 
sitic mode of a mesa beam cavity is only a factor ~ 2.5 
closer to the ground state frequency than the nearest par- 
asitic mode of a Gaussian-beam arm cavity. Therefore, 
as discussed in MBI Sec. IV C, we crudely expect the 
mesa-beam cavity to be only marginally more sensitive 
to perturbations than a Gaussian-beam arm cavity. 



The discussion of the previous section gives us good 
reason to believe that mesa-beam and gaussian-beam 
arm cavities will display the same sensitivity (i.e. within 
a factor ~ 2.5) to mirror perturbations. Here, we test 
this proposition when the ETM is displaced through a 
distance s = sx - or, more explicitly, when mirror 2 is 
perturbed by 5h2 



Shdisp, given by 



Adh2 
[s ■ r)- — 
r ar 

1 |s X dh2 

2 r-^ dr 



1 {s-r)^ d^h2 

2 r2 Ik^ 



^ ^ 1 dh2 ^ 1 {xs)"^ d^/i2 
r dr 2 r^ dr^ 



(56) 



+ 



1 (ys) dh2 

2 r^ dr 



Section II V El summarizes the step-by-step process we 
follow to explore the influence of perturbations (here, dis- 
placement). However, because we expect — and our calcu- 
lations below confirm — that whatever the precise mirror 
shapes, the properties of the interferometer will depend 
only weakly on displacement,^^ we do not complete all 
the steps that procedure includes (e.g., we do not com- 
pute the change in thermoelastic noise or cavity gain with 
displacement). 



^'^ The natural length parameter for the problem is fe = \/ XL/2tt = 
2.6 cm, the diffraction length. When we perform perturbation 
theory, we find results which vary in powers of s/b. Since the 
LIGO control system will control displacements to much smaller 
than 1 cm, displacements have a relatively small effect on the 
LIGO interferometer. 
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1. Displacement of fiducial spherical mirrors 

To evaluate the perturbation expansion [Eqs. (|32|1 and 
(|34|l ] for a cavity bounded by two identical fiducial spher- 
ical mirrors subjected to a displacement of its ETM 
through a distance s [Eq. JSHJ], we use analytic tech- 
niques specialized to spherical mirrors (e.g., Hermite- 
Gauss basis functions; cf. Appendix If)| . After some al- 
gebra (described in detail in Appendix IF|| . we find the 
ground state of the perturbed cavity, to first order, to be 
[cf. MBI Eq. (4.12)] 



sph 



io)+cr' 1(1,0)) 



V2(l + .g)3/4 
= 0.008(s/l mm) 



{s/b) 



(57a) 



(57b) 



[For clarity and for consistency with other work, we 
have absorbed a phase into the definition of the |(1,0)) 
Hermite-Gauss state.] Here, g =1 — L/Rc = 0.952 is the 
value for the two fiducial spherical mirrors. 



2. Displacement of fiducial Mexican-hat mirrors 

Similarly, if we repeat the above calculation for 
Mexican-hat mirrors using the methods summarized in 
Sec. IIV El we find the ground state of the cavity changes 
to [cf. MBI Eq. (4.12)] 



Wi) 
MH 



CI 



o)+crK 

0.263 {s/b) = 



)+0{<f) (58a) 
(58b) 

0.010 (s/1 mm) . (58c) 



Note that unlike the Gaussian-beam case, \'Wi) is a unit- 
norm superposition of eigenmodes of the unperturbed 
(mesa-beam) cavity, rather than an eigenmode of that 
cavity itself. 



3. Power in parasitic modes 

Perturbations cause changes in the resonant ground 
state. Equivalently, perturbations couple the ground 
state to the parasitic modes, causing power to leak from 
the resonant ground state of the unperturbed arm cavity 
into these other arm cavity eigenmodes. 

In the case of displacement, the power in the parasitic 
modes is easily distinguished from the carrier light by 
symmetry: while the carrier light is axisymmetric, the 
lowest-order changes in state are dipolar (i.e. [V'^^-') « 
1(1,0)) for Gaussian beams). The fraction of the total 
arm cavity light power in the dipolar parasitic modes is 
[cf. MBI Eq. (4.13)] 



2 f 100(s/1.3mm)^ppm sph, 
1 = Ci - I ioo(s/1.0mm)2ppm MH. 



(59) 



[In deriving the above expresison (|59|l . we assume that 
for each displacement distance s, the input beam has 
been optimally tuned to match the arm cavity's per- 
turbed ground state eigenfunction. In a sense, this ex- 
pression measures the intrinsic fraction of power sent 
into parasitic modes by the applied perturbation (here, 
displacement of one ETM). 

In practice, the arm cavity is driven by some fixed 
driving beam. As the mirror is displaced, the fraction 
|7o(s)|2 of input laser power that enters the arm cav- 
ity will change [cf. Eq. (|37|l ]. All this laser light enters 
the arm cavity's perturbed ground state. Of the power 
that enters the arm cavity, only a fraction « Pi enters 
parasitic modes. Therefore, relative to the total power 
beamed towards the arm cavity, only a fraction |7o|^Pi 
enters parasitic modes.] 



4. Power out the dark port 

If one arm cavity's ETM is displaced, then the light 
leaving the two arm cavities will not interfere destruc- 
tively at the dark port. The precise amount of power 
Pdp out the dark port depends on the driving beam [cf. 
Eq. H40|) . which depends on 70, 71, and 72]. If the in- 
terferometer is driven by the optimal Gaussian beam (cf. 
Sec, nil A^ . then the power out the dark port is approx- 
imately entirely in a dipolar mode, with net power [i.e. 
Eq. igOI) to lowest order, with 71 = 0; cf. MBI Eq. (4.14)] 



P. 



total 



dp 



( 90(s/1.3mm)2ppm sph, 
\90(s/1.0mm)2ppm MH. 



(60) 



C. Effect of tilt on the resonant eigenstate of an 
arm cavity bounded by spherical and Mexican-hat 



mirrors 



Because both the spherical and Mexican-hat mirrors 
planned for advanced LIGO are very flat, the planned 
advanced LIGO arm cavities will necessarily be fairly sen- 
sitive to tilt. 2* This section demonstrates that, though 
Mexican-hat mirrors have a very flat central region — 
much more so than their spherical counterpart (cf. Fig. 
|2Jl — an arm cavity bounded by Mexican-hat mirrors will 
only be somewhat (i.e. a factor ~ few) more sensitive to 
tilt. 

If the ETM of an arm cavity is tilted through an an- 
gle 0y about its y axis, the mirror surface is effectively 



For example, to order of magnitude, a tilt angle 9 ^ \J \/2ttL = 
6 X 10~^ should produce a perturbation of order unity in the 
optical state of the cavity. Because an order-unity change in state 
implies fairly dramatic change in the interferometer, tilts much 
smaller still (of order fewxlO"*) can cause serious difficulty with 
the interferometer. 
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perturbed by Sh2 = Shtnt- 

Shtm = Oyx + Oiel) 



(61) 



1. Tilt-induced changes in the arm cavity ground state 

When this perturbation is inserted into the perturba- 
tion expansion H32I) and the terms in that expansion are 
evaluated in the case of spherical mirrors (using special 
properties of spherical mirrors and Hermite-Gauss basis 
states; cf. Appendix IF|| . we find first- and second-order 
corrections to the state [i.e terms in Eqs. (|32() and (|34|l ] 
to be given by [cf. MBI Eqs. (4.3) and (4.4)] 



iV-i) 

sph 



1(1,0)) , 

s 

"2 



^^^P^|(2,0)) 



1 



-(0y X VkL) 

V2(l - .g2)3/4 ^ y ' 



O.OO64(0j,/lO"^) 



sph 



^/2(l-52)l/2(l_^) 
= 0.00046(6iy/10"^)2 . 



{9y X Vfcl)^ 



(62a) 
(62b) 

(62c) 

(62d) 



Here g = 0.952 is the g-value for the fiducial cavity, the 
states \{m,n)) denote states in the Hermite-Gauss basis, 
and these states have been adjusted in phase to make ai 
and a2 real. 

If this same expansion is evaluated using Mexican- 
hat mirrors via the numerical procedure outlined in Sec. 
El we find first- and second-order corrections to the 



H 

state [i.e terms in Eq. H34|l ]. We will not provide an ex- 
plicit form for these states here. These corrections have 
norms given by the Mexican-hat analogues of Eq. 1)62(1 
for ai and a2 [cf. MBI Eq. (4.5)]: 



,MH _ 



"2 



IIV-'^'ll = 14.78 {OyVkL) 

0.0227(0/10"*) 

IIV-'^^II = 74.97 (OyVkLf 

omoi8{ey/io-^f . 



(63a) 



(63b) 



We can use these expressions to define normalized repre- 
sentations of the first- and second-order corrections: 



Ml 
U2 



(63c) 
(63d) 



We display contour maps of these normalized admixtures 
of modes in Fig. 5 of MBI 



Gaussians; for mesa beams). Therefore, interpreting 
this change as an excitation of dipolar parasitic modes, 
the fractional power in the dipolar parasitic modes is [cf. 
MBI Eq. (4.6)] 



Pi 



al 



0.0005(6ly/3.5 X 10-8)2 
0.0005(61^/1.0 X 10-8)2 



sph, 
MH, 



(64) 



when the ETM is tilted through an angle 9y 



3. Tilt-induced changes in the diffraction losses of the 
ground state of the resonant arm cavity (Mexican-hat only) 

The diffraction losses associated with the ground 
state — which we approximate by the clipping approxima- 
tion losses £i and £2 [Eq- l|16|l ] — also change when the 
beam state changes. As described in Sec. lIVCTTl we find 
the perturbed value for, say, Ci merely by expanding the 
expression for Ci, obtaining the general expansion Eq. 
(El- 

If we evaluate Eq. H36|l for the case of a cavity bounded 
by two Mexican hat mirrors with one mirror (the ETM, 
i.e. mirror 2) tilted through an angle 9y, we find we can 
rewrite Eq. H36|l for the losses at mirror 1 (the ITM) in 
terms of the expansion mentioned above (i.e. in terms of 
ai, a2, ui, and U2): 



a2'C^' 



£'1 = Ci+al 

w 18 ppm [1 + 0.0025(6lj,/10-8)2] 

[cf. MBI Eq. (4.7)] where 



(65) 



Ci = 

,(2) _ 

-A ~ 

,(2) _ 



{uo\Oi\uo) = 18 ppm , 
(ui|Oi|ui) = 96 ppm , 
2Re(Mo|Oi|M2) = 29 ppm 



(66) 



and where ui and U2 are defined by Eq. I|63|) . 

Using a similar technique we can also evaluate the 
change in the clipping approximation diffraction losses at 
mirror 2. [As discussed Sec. lIVD2l we for simplicity as- 
sume the diffraction losses at both mirrors remained the 
same as the mirror tilted, or Ci = £2- Since the diffrac- 
tion losses influence physical quantities like the cavity 
gain fairly little, we require only a rough estimate of their 
sensitivity to perturbations.] 



D. Effect of tilt on arm cavities and 
interferometers using Mexican-hat mirrors but 
driven by Gaussian beams 



2. Parasitic mode power excited by tilt 

The largest correction to the perturbed arm cavity 
ground state is a dipolar perturbation (i.e. |(1,0)) for 



Tilt of one ETM causes changes in several important 
properties of an interferometer that uses mesa beams 
(and is driven by the optimal Gaussian): (i) the amount 
of thermoelastic noise present in the dark port signal; 
(ii) the amount of power present in the perturbed arm 
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cavity (i.e. the arm cavity gain); and (iii) the amount of 
input hght power leaving the dark port. In this section 
we compute the changes in these three quantities by ap- 
plying the general techniques presented in Sec. IIV CI to 
the case of tilt. 



1. Tilt and thermoelastic noise 

To lowest order in Oy/{h/L), tilt has no effect on ther- 
moelastic noise: at first-order, tilt only excites odd-parity 
modes, which produce nonaxisymmetric intensity pertur- 
bations, and nonaxisymmetric intensity perturbations do 
not contribute to lowest-order changes in the thermoelas- 
tic noise integral (cf. Sections IIV C 41 and IIV D T)) . 



2. Tilt and cavity gain, when driven by a Gaussian beam 

Equation (|38|l provides a general expression for the 
cavity gain for a perturbed cavity; this expression de- 
pends on (i) the amount of light power 7^ entering the 
perturbed arm cavity [Eq. %\7\ . evaluated using the per- 
turbed mesa-beam state from Eq. Hfi3|l ] . which we can 
express using the perturbation parameters 71 and 72, 



71 = 72 = -0.070 + 0.013i 



(67) 



and (ii) on the diffraction losses of the resonant state of 
the tilted arm cavity [Eq. 1)65(1 above]. Combining these 
expressions, we find the arm cavity power, relative to the 
input beam power, to be given by [cf. MBI Eq. (4.8)] 



anal 



l-o? + ^(72 +72)- 2^^ 

70 1 - r/ 



= 737 [1 - 5.5 X 10"^(6'j,/10"*)2] 



(68) 



Tilt and dark port power, when driven by a Gaussian 
beam 



Finally, when one ETM mirror in a mesa-beam inter- 
ferometer is tilted, the beamsplitter sends light to the 
dark port. Equation l|40|) provides a general expression 
for the light Udp leaving the interferometer through the 
dark port; for tilt, this expression evaluates to 

Mdp = \/2 [(-70Q:? + 720:2) \uo) + ai7o jwi) + a27o |m2>] 

where we use (i) the definitions ai = | \ \ and similarly 
[cf. Eq. (|^ ]. and (ii) 71 — 0, to simphfy the general 
equation 



The corresponding fraction of the interferometer's 
power that exits the interferometer through the dark port 
in the fundamental mode Uo and in the parasitic modes 
Ml and U2 is [cf. MBI Eq. (4.9); note the MBI expression 
will be larger by a factor 4]"^° 



-Pdp,0 

-Pdp,i 

-Pdp,2 

E. 



loal 

MH\2 



7o ("l ) 



72a2|2 ~ 0.256(V10"**)'^ Pl(iOa) 
2 _ (70b) 



478(6lj,/10-'*)- ppm 



%{af''Y ~ O.O24(0j,/lO-«)^ ppm . (70c) 



Effect of mirror figure error on mesa-beam 
interferometers 



In this section, we explore the sensitivity of individual 
arm cavities and the overall interferometer to mirror fig- 
ure error, when the beam resonating in the arm cavities 
is a fiducial mesa beam. More specifically, in this section 
we (i) distort the ETM by a physically plausible amount 
(i.e. an amount estimated from actual mirror figure error 
measurements of LIGO mirrors); (ii) apply perturbation 
theory (cf. Sec. IIV E|) to deduce the change in resonant 
ground state of the arm cavities; and then (iii) use the 
resulting modified beam state to deduce how the power 
spectrum of thermoelastic noise will change due to the 
perturbed beam cross section. 

MBI i uses the computations performed in this sec- 
tion to place constraints on the accuracy of machining 
required of the mirrors used in a mesa-beam interferom- 
eter; cf. MBI Sections IV F and IV G. 



1. GLB's worst-case figure error 

GariLynn Billingsley has provided us with a map of a 
worst-case figure error, 5z^c(x,y) [height error as func- 
tion of Cartesian coordinates in the transverse plane], 
produced by current technologies. Her map is based on 
the measured deviation of a LIGO-I beam-splitter sub- 
strate from flatness. The measured substrate had di- 
ameter 25 cm; she stretched its deviation from flatness 
(its "figure map") to the baseline advanced LIGO mir- 
ror diameter of 35.4 cm, fit Zernike polynomials to the 
stretched map, and smoothed the map by keeping only 
the lowest 36 Zernikes. MBI Fig. 6 provides a contour 
diagram of the resulting figure map (figure "error"). In 
the central region (innermost 10 cm in radius), the peak 
to valley error Az is about 30 nm, while in the outer 
region (10 cm to 16 cm in radius), it is about 110 nm. 

Billingsley (private communications) thinks it likely 
that in the central region (which dominates our consid- 
erations), peak-to- valley errors of Az ~ 5 nm (about 1/5 



When computing the diffraction losses for a tilted cavity, we 
assume the same diffraction losses associated with a bounce off 
mirror 1 and mirror 2. 



The MBI expression estimates the effect when all four mirrors 
are tilted about uncorrelated axes; it therefore is larger than our 
result, which describes the effect of tilting only one mirror. 
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weight fiiTiction ,[cf. Mm Kg, (4-1 9)]: 



e,W {r/b - 3.7) 5z^, 

[1-W {r/b - 3.7)] Sz^ 




(72) 



> 1, W{x) = X if X e [0,1], and 
is tapering function permits us 
n in radius one level of figure 
the outer ~ 6 cm a different 



-■without compensating tilt 

lionted according to Sec. 
I the change in the mesa- 
2 is perturbed by GLB's 
Id by a factor e [E g. ((7T|) ]. 
])h from Sec. IV CI to eval- 
cond-order corrections to 
HEgs. (4.16) and (4.17)] 



•2) 

= 0.167e 
0.123e2 



15 -10 







(73a) 
(73b) 
(73c) 
(73d) 

sciiifi^taet^^Bl^st be fairly small (~ 0.2) for the 
secojid'c)j;d:ffl>^^BJ^pion to be significantly smaller than 
efk^^^j^PErrection (i.e. for the series to converge, 
aiidr-'^^B^^wlaation theory to -be useful).] The series 
p'can be used, as in the case of displacement 
and tilt, to determine physical guantities such as 



The powlrQn parasi in the arm cavity: 

MBI Fig. 7 shows a map of the power distribution 
|^(i)(7.)|2 _ i^^y^^-^-jp jj-^ -j-j-^g most-strongly-excited 

parasitic mode combination \vi). The net power in 
parasitic modes (as a fraction of total power in the 
arm cavity) is well approximated by the integral of 
this guantity [cf. MBI Eg. (4.17)]: 

Pi = ^ 0.0011(e/0.2)2 . (74) 

The power leaving the dark port: Eguation l|40|) 
tells us that the fraction of the interferometer's 
overall laser light input power going out the dark 
port is approximately [i.e. Eg. (|40|) to lowest order; 
cf. MBI Eg. (4.18)] 



X, cm 



FIG. 5: Contour diagram of GariLynn Billingsley's worst-case 
mirror figure error [height Szv,c] in nanometers, cf. MBI Fig. 
6. 



as large as in MBI Fig. 6) may be achievable; Jean Marie 
Mackowski believes even smaller errors can be obtained 
with coating methods (cf. MBI Sec. IV F 1). Accordingly, 
in the analyses described below we shall use Billingsley's 
map, scaled down in height by a factor e [cf. MBI Eg. 
(4.15)]: 



6z = e6z„c{x,y) , 



(71) 



Pdp « (l7o|' + |7i|')|/3i|' 

~ 72[/3i[2 =.0.0010(6/0.2)2 



(75) 



and we shall use e = 0.2 and Az = 6 nm as our fiducial 
values for e and Az. 

Also, Billingsley thinks it likely that the outer regions 
of the mirror will be of significantly lesser guality than 
the inner regions. To study the sensitivity of the interfer- 
ometer and thermoelastic noise to errors in the exterior, 
we divided GLB's perturbations into two regions by a 



where we neglect 71 and compute only the lowest 
order term in the dark port power. 



Unlike displacement and tilt, for mirror figure error 71 is gener- 
ically nonzero, since the perturbation admits an axisymmetric 
part. 
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Naturally, the interferometer is more sensitive to er- 
rors in the inner ~ 10 cm, where the beam power is large, 
than to the outer 6 cm, where the beam power is small. 
To investigate this effect, we considered height perturba- 
tions of the form discussed in Eq. (|72|l : to characterize 
the sensitivity of the beam shape to these defects, we ex- 
amined the norm = [Eq. H73|) ]. which is neces- 
sarily a quadratic form in Sc and Eq [cf. MBI Eq. (4.20)]. 
Using the same methods as for Eq. (|73|l . we found the 
norm Pf to be 

Pf = OmSel + 1.06 X 10-^eo{eo - Sc) (76) 
« 0.0011(ec/0.2)2 -I- 4.3 x IQ-^ {e o/0. 2) {eo - ec)/0.2 . 

If P\ is used as a characteristic example of the sensitivity 
of physical quantities (e.g., the power out the dark port) 
to height perturbations, then a mesa-beam interferom- 
eter is around 5 times more sensitive to mirror figure 
errors in the interior region of each test-mass mirror (i.e. 
its inner ^ 10 cm) than it is to perturbations outside 
that region. ■^^ 



3. Mirror figure errors, with compensating tilt 

The tilt control system automatically and dynamically 
reorients the mirrors in response to what it interprets as 
tilt. Specifically, the mirror tilt control system (i) mea- 
sures signals containing information about beam asym- 
metry^ such as the output of a quadrant photodiodc; (ii) 
computes the mirror tilt that would generate these asym- 
metries using perturbation theory expansions [i.e. Eq. 
(|32|l . to first order]; and then (iii) tilts all four LIGO 
mirrors to eliminate the apparent tilts that the system 
computed in the previous step. Therefore, since mir- 
ror figure error also produces beam asymmetries, the tilt 
control system of the interferometer will act to partially 
compensate (the dipolar part of) the static mirror defect. 

The precise quantity the tilt control system measures 
to deduce the tilt angle is not important: different ap- 
proaches to tilt control interpret the optical state of our 
arm cavity [i.e. interpret the state defined by the 
corrections in Eq. H73|) ] in a manner fairly independent 
of the method used (i.e. the compensating tilt we calcu- 
late depends little on the method we use to calculate it). 



In this section, we assume the tilt control system acts to 
minimize the dipolar component of the arm cavity beam 
power. 

Therefore, wc find the tilt-compensated state by (i) 
adding together the results of a tilt perturbation [i.e. a 
generalization of Eq. H62|l that depends on two tilt an- 
gles. Ox and 9y\ and the mirror figure error perturbation 
[Eq. H73|l . which depends on e]; (ii) evaluating the norm 
I IV'^^^ 1 1 of the first-order perturbation; and then (iii) find- 
ing the tilt angles 9^ and 9y which minimize that norm. 
The resulting optimal tilt angle is [cf. MBI Sec. IV H 3] 

9x = +0.98 X 10"^(e/0.2) , (77a) 
9y ^ +0.69 X 10"^(e/0.2) , (77b) 

e = ^91 + 91^ 1.2 X 10-^(e/0.2) (77c) 

(in the limit of small e, so linear theory applies). MBI 
Figure 8 shows the height of the surface of the mirror 
after the compensating tilt is applied. 

After the tilt is apphed, the first-order correction to 
the resonant mesa-beam state is a combination of tilt and 
mirror-figure perturbations. For brevity, we denote the 
net first-order correction to the mesa-beam state when 
tilt is applied by 

V-^'^) = Pi,cM , (78a) 
= ||V'(i)||0.02(£/0.2). (78b) 

The square of the norm (i.e., Pf as usual, is the frac- 
tion of arm cavity power which is in parasitic modes [cf. 
MBI (4.23)]: 

Pf™ = |/3i ,|2 = 0.0004(e/0.2)2 . (79) 

MBI Figures 9 shows a map of the power in parasitic 
modes, \Pi,cVi^c{t)\'^ . Also, the norm of this first-order 
correction provides an estimate of the power leaving the 
dark port of this interferometer (i.e. an interferometer 
with one tilted, defective mirror) when the interferome- 
ter is driven with Gaussian beams [i.e. Eq. (|40l to lowest 
order; cf. MBI {4.24)]-^ 

Pdp « (|7o|' + l7i|')l/3i,c|' 

-fllPif ^ 0.00038{e/0.2f (80) 



^■^ In our perturbative computations, the mirrors were treated as in- 
finite. The infinite-mirror model proves sufficiently accurate for 
perturbations in the center 10cm of the mirror. Unfortunately, 
for mirror defects in the outer 6cm, the resulting perturbed state 
has diffraction losses roughly comparable to the effect of the per- 
turbation. Therefore, our perturbative simulations of mirror de- 
fects in the outer ~ 6cm in general, and the coefficients of to in 
Eq. in particular, are not reliable. By contrast, in her sim- 
ulations Erika D'Ambrosio used realistic finite mirrors 0. As 
discussed in MBI [cf. MBI Eq. (4.20)], she finds the inner region 
to be about 3.5 times more sensitive to mirror perturbations than 
the outer region. 



4-. Influence of mirror figure errors on thermoelastic noise 

Because the ETM's figure error distorts the beam res- 
onating in the arm cavity, the thermoelastic noise pro- 
duced by each mirror bounding that arm cavity changes 
by some small amount. Given the change in state de- 
duced above, we know how the beam profile at the ITM 



As in the untilted case, we neglect the (nonzero) term |7ip. 
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changes. We can therefore evaluate, using the discussion 
of Sec. II V 01 the change in thermoelastic noise associated 
with the beam reflecting off the ITM's face. 

To be explicit, to compute the first-order effects of a 
perturbation 5P to the thermoelastic noise integral /i of 
mirror 1 in this fiducial case, we first select the axisym- 
metric portion 5Po{r) = / dipP{r) /2n of 5P; we then lin- 
earize the analytic expressions derived by Liu and Thorne 
for cylindrical mirrors (cf. Sec. IV A"2|l . using 5Po as the 
magnitude of the perturbation'^'*; and finally we extract 
from this linearization the first-order change in Ii . 

When we apply this technique to the fiducial beam re- 
flecting off a cylindrical mirror of radius 16 cm and thick- 
ness 13 cm which is deformed by GariLynn Billingsley's 
mirror distortion [Eq. H71|l ]. we find the thermoelastic 
noise integral for that mirror changes to 

/i/(2 X 10^28) ^ Q [1 0.035(£/0.2)] (81) 

where /i denotes the total thermoelastic noise integral 
for mirror 1 when mirror 2 is deformed by our scaled 
height perturbation l(7T1l . [Roughly speaking, we expect 
the beam and hence the thermoelastic noise to change 
in a qualitatively similar fashion at mirror 2 (cf. Sec. 
llVD2t ]. 

MBI Section IV G applies this result to deduce how 
sensitive the power spectrum of thermoelastic noise is to 
uncorrelated mirror figure errors on all four mirrors [cf. 
MBI Eq. (4.27)]. 

VIII. CONCLUSIONS 

In this paper, we have described both the theory and 
practice needed to obtain the results summarized in MBI. 
Specifically, we have developed analytic and numerical 
tools to evaluate the following: 

• Thermoelastic noise integrals: We developed prac- 
tical techniques for finding the thermoelastic noise 
for nonstandard optical systems (i.e. noncylindrical 
finite mirrors and unusual beam shapes) [Sections 
III Bl and IV Aj . We tabulated the thermoelastic 
noise, relative to the current advanced LIGO base- 
line, for many alternative mirror and beam con- 
figurations [Tables HI -HVl in Section |^. We found 
many configurations with lower thermoelastic noise 
than the baseline advanced LIGO configuration. 

• Eigenmodes of an arm cavity hounded by Mexican- 
hat mirrors: We also wrote numerical code to find 
the eigenmodes of an optical cavity [Sections III Gl 
and IV B| . We computed and tabulated many of the 



eigenmodes of an arm cavity bounded by Mexican- 
hat mirrors (none of which appear explicitly in this 
paper). 

• Second- order optical perturbation theory: Finally, 
we developed expressions for second-order optical 
perturbation theory [Sec. llVj. We applied pertur- 
bation theory extensively to study the sensitivity of 
mesa-beam interferometers to perturbations (i.e. to 
mirror figure error, mirror tilt, and mirror displace- 
ment). 

The results found here are used in MBI to conclude that 
mesa-beams interferometer designs offer clear advantages 
over the baseline advanced LIGO design, without being 
substantially more sensitive to mirror figure error, tilt, or 
displacement perturbations. 

In MBI U, in Erika D'Ambrosio's paper 6], and this 
paper, we and our collaborators have only taken the first 
steps towards the design of a practical mesa-beam ad- 
vanced LIGO proposal. For example, more perturbative 
calculations — this time, applied to the final design, rather 
than to a fiducial case — are needed, for the design of the 
control system (i.e. so the relationship between light on 
various photodiodes and the correcting tilt applied to 
the LIGO mirrors can be established). Further, in this 
paper we have only begun to explore the space of all 
possible mirrors. For simplicity, we chose to fix the mir- 
ror mass to 40 kg. In practice, however, with sapphire, 
mirror designs are limited by fabrication limits (i.e. the 
radius of the mirror is limited by the radius of the sap- 
phire boule one can grow) rather than by weight limits. 
Therefore, before a final design is chosen, more mirror de- 
signs (including cylinders with m > 40 kg) should be ex- 
amined. Finally, in this paper, our Mexican-hat designs 
were limited to symmetric cavities (i.e. using identical 
mirrors). In practice, asymmetric designs offer the pos- 
sibility of lower thermoelastic noise and greater practical 
convenience'^^. Further work is necessary on the design, 
construction, and operation of asymmetric cavities with 
mexican-hat-like mirrors. 
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APPENDIX A: SAPPHIRE MATERIAL 
PARAMETERS: NOTATION AND VALUES 

While sapphire is not an isotropic material, it can be 
reasonably approximated by isotropic elastic and ther- 
modynamic properties. In this paper, we therefore treat 
sapphire as an isotropic medium, with the following spe- 
cific values for physical parameters:^^ 



E 


= 4 X lO^^erg/cm^ 


(Al) 


P 


= 4.0 g/cm'^ 


(A2) 


a 


0.28 


(A3) 


a 


= 5.5 X IQ-^R-^ 


(A4) 


K 


= 3.3 X lO^ergcm-^K^^s^^ 


(A5) 


C'v 


= 7.7 X 10*^ergg-iK"i 


(A6) 



APPENDIX B: GENERAL DESCRIPTION OF 
MESA BEAMS 

In Section III Al we only briefly describe the form of 
mesa beams evaluated at the surfaces of an arm cavity's 
two mirrors. In this appendix, we provide a much more 
thorough treatment of mesa beams. 



1. Constructing mesa beams 

Mesa beams are constructed by averaging some gaus- 
sian beams (with waist size Ao and waist location z^) 
over some disc of size D and then normalizing the result. 



Specifically, we construct mesa beams that propagate to- 
wards positive z by first averaging a (non-normalized) 
gaussian beam over a disc of radius D: 

U{f,z,D) = 6*'=^ 

X / d?r'Ug{r ~ ? ,z;z,Zw)Bl&) 

Jr'<D 



Ug{f,z;z,Zw) = cxp 



1 - i 



X e 



-i tail" "'" (js/ z) 



i2 



(Bib) 



[where Aq is a function of z, given by z = kA^, and 
where Zw is the location of the beam waist]. The re- 
sult of the average is necessarily axisymmetric. Since 
C/(r, z, D) is axisymmetric (because of the symmetric in- 
tegral which defines it), we may normalize the result with 
a one-dimensional integral: 

/•oo 

N^(D,z) = / \U{r,z,D)\^2Trrdr (Blc) 
Jo 

u^,s.{r,z,D) = U{r,z,D)/N{D,z)e'''' (Bid) 

To construct the associated Mexican-hat mirror which 
will reflect this mesa beam back into itself (propagating 
in the opposite direction) at location Zr, we require the 
Mexican-hat mirror height function ft-MH to be continuous 
and to satisfy 



i(r, Zr, D) exp(— 2z/c/i 



MH) 



,(r,z„i?)* (B2) 



for ~ u{r, Zr, D) the values of a mesa beam at the mir- 
ror plane z = z,.. This expression is equivalent to MBI 
Eq. (2.13) [cf. Eq. JB)) in this paper]. These requirement 
uniquely specify the Mexican- hat mirror shape /imh- 



a. Canonical mesa beams: mesa beams for symmetric 
cavities 

The mesa beams presented in MBI and discussed 
henceforth in this paper — denoted canonical mesa beams 
— are assumed to have particularly special form. First, 
it is assumed that the cavity confining the mesa beam is 
symmetric, with mirror 1 placed at z = —L/2 and mirror 
2 placed at z = i/2; also it is assumed that the waist lo- 
cation, Zu,, is placed precisely between them (i.e z^, = 0). 
Second, it is assumed that the gaussian used to construct 
the mesa beam is the minimal- diffraction gaussian for the 
symmetric cavity. The minimal diffraction gaussian has 
z = L/2 and Zw = 0, implying that the beam width is 
b = \/ L/k at the mirror surfaces and 6/\/2 at the beam 
waist at z = 0. 



2. Canonical mesa beams at the mirror surfaces 

These numbers are taken directly from the advanced LIGO 

project book Q, and the advanced LIGO summary web page In this paper and in MBI, we Only need to know the 
13. form of canonical mesa beams at the mirror surfaces. In 
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this case, since the beam is canonical (i.e. has symmetric 
cavities and minimal diffraction gaussians, so z^, = and 
z = implying Aq — h/y/2) and since we evaluate it 
at the mirror surfaces z = we find that Ea. (|Blb|) 

simplifies to 



5, mm 



(f,z = Ti/2) 



exp 



[l±^] 



262 



(B3) 



Inserting this expression into the definition (|BI|I of gen- 
eral mesa beams, we find that the mesa beam form at 
mirror I (i.e. the mirror at z = — L/2) is given in terms 
of the construction described in Eq. (jSJ. [For simplicity, 
in Eq. Q we omit the overall phase factor exp(i7r/4).] 



APPENDIX C: CONVERTING THE 
FLUCTUATION-DISSIPATION MODEL 
PROBLEM TO A STATIC MODEL PROBLEM AT 
LOW FREQUENCIES 



In Sec. Ill Bl we claim that at low frequencies, the elas- 
tic response of a cylinder to an oscillating pressure profile 
can be reconstructed to a good approximation using the 
static response of a cylinder to the same pressure profile, 
in an appropriate accelerating frame. 

Briefly speaking, a relationship between the dynamic 
and static problems exists because, when the oscillations 
are sufficiently slow, the effect of dynamical terms in the 
elastic equations of motions can be neglected. These 
physical considerations have been discussed at greater 
length elsewhere, cf. Sec. II of LT [l^ . 

In this appendix, we demonstrate in more technical de- 
tail precisely how to establish the desired (approximate) 
relationship. We assume only basic familiarity with elas- 
ticity, on the level of Blandford and Thorne 



; y w i 

El. 



1. General quasistatic approach for the elastic 
response to an oscillating surface stress 

We wish to solve the elasticity equations, which can be 
expressed as 



pdfya + y'Tab = , 



(CI) 



for a test-mass mirror subject to an oscillating surface 
stress on its inner face (i.e. on the surface z = 0, where 
the optic axis of the advanced LIGO arm cavity is the z 
axis) , 



Ta,{t, r,z^O) = Ta,{r, z = 0) cos{ujt) 



(C2) 



and otherwise subject to no other stresses. [Here, Tab is 
given in terms of ya by Eq. (|f Ob|) .] 



a. Step 1: Express the problem in the accelerated frame 

We can better understand the response of the mirror 
substrate if we go into a frame comoving with the center 



of mass of the test mass mirror. The test mass experi- 
ences a net force FaCOs{ujt), determined by the surface 
stress: 



cos(wt) / <frTaz{r, z = 0) . (C3) 



In response to this net force, the mirror center of mass 
R cm accelerates: 



R cm 



Fa 



cos(aji) 



(C4) 



Therefore, to go to the comoving frame, we perform 
the following transformation: (i) define the comoving dis- 
placement field y'^ = ya — Rem- 



y'a = ya + COS(WO 



(C5) 



(ii) define the comoving stress-energy tensor as T^f^ — 
Tabiy')- In terms of these two new quantities, the elas- 
ticity equations in the accelerated (primed) frame are 



M 



(C6) 



The right hand side is simply the inertial force associated 
with working in an accelerated frame. 

The boundary conditions remain unchanged. 



h. Step 2: Factor out all sinusoidal dependence 

In the accelerated frame, we now assume all quanti- 
ties oscillate sinusoidally in response to the sinusoidally- 
oscillating pressure profile and inertial force: 

y'a = Va cos(wi) , 

T'ab = fabCOs{ujt) . 

Substituting these expressions into Eq. (|C6|I . we find we 
can reconstruct a solution to the dynamic solution by 
solving the static partial differential equation 



2- I V7&rfi Fa 
-PUJ ya+y Tab = -Pj^ 



(C7) 



subject to force-free boundary conditions on all surfaces 
except the top surface, which is subjected to a constant 
pressure profile Taz{r, z = 0). 



c. Step 3: Approximate the problem as static in the 
accelerated frame 

Finally, at sufficiently low frequencies (i.e. frequencies 
so low that sound crosses the cylinder many times within 
one period, as is the case for advanced LIGO; cf. notes 
El and El, the first term in the accelerated- frame elastic- 
ity equations [i.e. the term oc y^uj"^ in Eq. (|C6I) ] can be 
neglected (cf., e.g., LT Sec. II). The remaining problem 
[i.e. Eq. (fTn|l ] can be interpreted precisely as equations 
for the static elastic response of a solid to an imposed 
pressure profile. 
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d. Summary 

To summarize, then, we find an approximate solution 
for the elastic response of a solid to an imposed surface 
stress Taz at z = by 



Ua = cos(a;t) 



^ Fa 



Va 



(C8a) 



where ya is obtained as a solution to the static elastic 
equations in an accelerated frame 



not change 6 (i.e. 6 = 9') and since (cos^(a;i)) = 1/2, 
we conclude that 

where the approximation neglects terms which are small 
in the quasistatic limit (cf. notes |21 and |2) . 



APPENDIX D: ON THE COMPLETENESS OF 
BASIS STATES FOR AN ARM CAVITY WITH 
TWO IDENTICAL, INFINITE MIRRORS 



(C8b) 



subject to the effective static surface stress Taz and bulk 
acceleration FajM given by 

Taz - Ta,COs(c^t) (C8c) 

Fa = - J d^rfazir,z = 0) (C8d) 

To summarize our conclusions on physical grounds, 
because the elastic response occurs much more rapidly 
than the surface pressure profile changes (i.e. because the 
sound crossing time is much shorter than the oscillation 
period of the imposed force) , we can effectively treat the 
elastic response as instantaneous. The elastic solid moves 
slowly through a sequence of static configurations. 

This quasistatic approximation, however, must be per- 
formed with care. If one neglects the dynamical terms 
entirely, as other authors working on this and related 
subjects have done then one finds elasticity equations 
without the bulk acceleration term. That equation is 
inconsistent with the static boundary conditions we im- 
pose. 



2. Using the quasistatic elastic solution to simplify 
our thermoelastic noise integral 

We can apply the quasistatic elastic solution we just 
developed [Eq. (|C8|l ] to find the thermoelastic integral 1^ 
[Eq. 0] associated with our specific elastic model prob- 
lem, where the surface stress imposed has (i) Taz = 
unless a ~ z and (ii) Tzz = — cos{ujt)FoP(r) [Eq. 
(|lUc|) ]. Since the accelerated-frame transformation does 



Because of the fluctuation-dissipation theorem, many other au- 
thors working on thermal and th erm oelastic noise (e.g., Cerdonio 
and Conti [271 : Liu and Thorne Il2l : ...) must solve a similar or 
identical elastic problem to deduce the effects of noise. And most 
choose to approach it using a similar quasistatic approximation, 
ignoring any dynamical effects. However, most make the qua- 
sistatic approximation before they go to the accelerated frame, 
rather than after. Only Liu and Thorne have correctly accounted 
for the effects of acceleration, though they add those effects in 
by hand later. 



In this appendix, we demonstrate that the eigenprob- 
lem for two identical infinite mirrors [Eq. (|13|l ] admits a 
complete set of orthogonal eigensolutions. In the text, 
we use the resulting eigensolutions as a basis for building 
perturbative expansions. 

Simply, the eigenequation admits a complete set of or- 
thonormal states because it is an eigenvalue problem for 
a unitary operator. Unitary operators always admit a 
complete set of orthogonal eigenvectors. The technical 
details for the infinite-dimensional (L^) proof we need 
here are far beyond the scope of this paper. However, 
the reader can understand the result by two simple ar- 
guments: 

Orthogonality of eigenvectors for unitary matricies I: 
Via hermetian intermediary: From any unitary operator 
we can define a hermetian operator H — ilnU (e.g., by 
a series expansion). The hermetian operator admits a 
complete set of orthogonal states. Hence so does U = 
exp{—iH). 

Orthogonality of eigenvectors for unitary matricies II: 
Directly: Alternatively, a familiar argument can be ap- 
plied to demonstrate that, if f7 is a unitary matrix, the 
eigenspaces associated with distinct eigenvalues are mu- 
tually orthogonal. 

Let Xa and \a) denote all the eigenvalues and (nor- 
malized) eigenvectors of a unitary operator U. Since 
U^U — 1, we know the eigenvalues are pure phase: 
{a\WU\a) = |Aa|2 = 1. Further, Since U-^ has \a) 
an eigenvector with eigenvalue I/Aq = A*, we can show 
(6| U = \b {b\ and thus 



= {Xa- Ah) {b\a) 



(Dl) 



Therefore, eigenvectors associated with distinct eigenval- 
ues are mutually orthogonal. 

Furthermore, in each degenerate subspace, we can di- 
agonalize to find an orthogonal basis. Therefore, for uni- 
tary matricies, all eigenvectors can be assumed orthonor- 
mal. 



APPENDIX E: SECOND-ORDER 
PERTURBATION THEORY FOR EIGENSTATES 
OF AN INDIVIDUAL ARM CAVITY 

In this appendix, we develop second-order perturba- 
tion expansions which relate changes in shape of one spe- 
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cific mirror of a symmetric cavity (the ETM) to changes 
in the resonant states in general (and to changes of the 
the ground state in particular) of that cavity. 

Since the precise form of the eigenequation depends 
on how one chooses to represent the state (i.e. at what 
plane, in what direction), many alternative and equiv- 
alent perturbation expansions can be derived'^^. The 
form presented in this appendix is that developed by 
O'Shaughnessy; Sergey Strigin and Sergey Vyatchanin 
performed their calculations using an independently- 
derived (but provably equivalent) approach. 

As in the text, in this section we make heavy use 
of standard quantum-mechanics operator notation for 
states and inner products, a notation described briefly 
in Section COKcf., e.g., |ili[l3)- 



1. Preliminaries: Setting up notation for the 
expansion 

We study the effect of changes in height of mirror 2 
on the solutions to the eigenequation (|13|l in the case 
of a symmetric cavity with infinite mirrors. For clarity 
and simplicity of notation, we redefine the eigenequation 
problem we perturb still further, into the following ex- 
pression: 



G 
T 

52 



= GTG\^) . 
= e-'''^GiG^ 
= e-'^^ = 1 4 
= 2kSh2 . 



ST 



(Ela) 
(Elb) 
(Elc) 
(Eld) 



ries, giving 



{Vo + em + e m + ■■) (V'o + eV'i + • • ■) (E2) 

^G{l + 5T)G (^o + eV-i + , 
and then matching orders on both sides; in this expres- 
sion, £ is a formal perturbation parameter, added to 
rescale the change in the reflection operator: 5T{e) = 
eST. When we work out order matching, we find (for 
P a projection orthogonal to the ground state Vo, also 
denoted |0)): 



V - ??o + (0 0) 



(E3a) 



STP- 



1 



-P5T 



\0) + —^P5f\Q) 



, 



ilo-G^' 



(E3b) 



1 



Vo-G^ 



P 



6T~(0 



ST 



Vo-G-^ 



PSf\0) 



where the operation O on an operator O is defined by 



O = GOG 



(E4) 



Here, 6h2 is the change in the (inward-pointing) height of 
the far mirror (mirror 2) and ST = T — 1. [The operator 
T, used only in this section, should not be confused with 
the truncation operators defined in Eq. Hll|) .] 



Perturbation theory expansion, expressed using 
operators 



We can furthermore substitute into the above expression 
the expansion 



6T ^ G{ -162 - T^^s ) G = ~iS2 ~ ^GSlG 



We construct a perturbation theory expansion by ex- 
panding both sides of the eigenequation Eq. ljEl|l in se- 

^ The relationship between these expressions need not be trans- 
parent. For example, one can represent the eigenequation using 
states represented "halfway" through a reflection of one of the 
mirrors (a choice which conveniently renders the resulting equa- 
tion always perfectly symmetric). The transformation between 



to give us the final form of the second-order expansion of 
the state: 

our representation and this one depends on the mirror height. 
Thus, when the mirror heights are perturbed, the relationship 
between these two representations involves a unitary transfor- 
mation that depends on the perturbation. 



J 



0)^^{0\G6iG\0)~{0 



S2P 



Vo-G-' 



Pdo 



Vo~G^' 



2 Vo - G2 



PGSlG\0) 



1 



Vo'G^ 
I 



P 



Vo-G^ 



PS2 |0) + 0(61) 



3. Perturbation theory expansion, expressed using rather through a finite collection of matrix elements of 
basis states 



As a practical matter, we compute the perturbation 
series expansion using not the operators themselves, but 
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the relevant operators relative to basis states. Therefore, 
we insert the identity operator, represented as a sum over 
all basis states (i.e. 1 = J2k 1^) (^1)' several points in 
the above expression. Since the basis states are eigenvec- 
tors of the propagation operator G and since the operator 
G is unitary, we know 



G\k) 
{k\G 



Xk\k) 
Xk {k\ 



(E6a) 
(E6b) 



Therefore, we conclude that the eigenvalue changes as 



r; - r;J l-z(0|<52|0) 



4(o|^i|o)-E 



fe#0 



%l(0|'^2|fc)| 
Vo - Vk 



(E7a) 
O (SI) 



while the state changes according to the expansion 



Vo - Vk 



(E7b) 



-?(H'52|0) + -^(fc|<52|0) (0\6,\0)-J2-^{k\S2\p){p\S2\0) 



Vo - Vk 



Vo - Vp 



APPENDIX F: PERTURBATION THEORY FOR 
CAVITIES BOUNDED BY TWO IDENTICAL 
SPHERICAL MIRRORS 

For cavities bounded by spherical mirrors — that is, 
mirrors with height function /11.2 = /2TZi^2 — the eigen- 
functions are known and of simple, tractable Hermite- 
Gauss form. We can therefore perform perturbation- 
theory calculations analytically. In this appendix, we 
describe these basis states and their application to per- 
turbation theory in greater detail. 



1. Background: Notation and definitions 

A spherical mirror of height h (r) is uniquely charac- 
terized by its radius of curvature TZ: h — /2TZ. We con- 
sider two identical such mirrors, placed symmetrically at 
zLL/2. We drive this cavity so it has kL/2'K ^ 1 wave- 
lengths between the mirrors (so kL is to an excellent 
approximation independent of mirror shape). 



a. Coordinate representation of eigenfunctions 

The eigenfunctions for a cavity bounded by two iden- 
tical spherical mirrors are known 15, IGI IitL flM Il9l | : 



Here Hm are Hermite polynomials, Ag = \/ z/k, N = 
m + n, g = 1 — L/TZ, 9 = cos^^ g, k = {qn + 0)/L for q a 
large integer, and 



A{zf ^ Al[l + {z/zf) 



2V 1-5 



(F2a) 
(F2b) 



In particular, at the mirror faces at z± = ±L/2, we 
have 



z±/z 



± 



l-.g 
1 + 5' 



(F3a) 



tan-^(z±/z) = ±-cos"ig = ±61/2 , (F3b) 



A{z±) 



fc (1 - .g2)i/4 



(F3c) 



6. Basis as eigensolutions of cavity 

After some algebra, one can verify these states iprn,n 
are solutions to the round-trip eigenequation Eq. (|13|l [or 
equivalcntly Eq. ^] with hi{r) = h2{r) = r^/2TZ, and 



\ — fP, — -ie{m+n+l) 

^^'mn V 'Iran — ^ 



(F4) 



1 



2™+"m!n! 
1 



H„ 



A{z) 



Hn 



A{z) 



^A{z 
X exp 



■ exp 



2A{z) 



i {N + l)tan~^ - -I- ikz 



(Fl) 



c. Establishing a quantum-mechanical correspondence 

For technical reasons not discussed further here, the 
eigenfunctions of a cavity bounded by spherical mir- 
rors correspond directly to the states of a 2-dimensional 
quantum-mechanical scalar nonrelativistic particle in a 
quadratic potential (i.e. a 2-d quantum simple harmonic 
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oscillator, or SHO). Like the SHO, the states are highly 
degenerate. If we index the states using their cartesian 
symmetry properties [i.e. express them as a Hermite- 
Gauss basis, as in Eq. ljFl|l ]. these quantum states of 
a SHO have the form |r7i,rt) = ^'m(a;)5'„(2/) for 



' ' ,HJ^)e-^y^-' (F5) 



for A some length scale. Therefore, if we evaluate the 
fields at the mirror surface z = —L/2 and define A = 
A{-L/2), we find 

= -L/2) - e'*"*™(x)*„(y) (F6a) 

where is given by 

= (N +1)9/2 + (independent of N) . (F6b) 

Since the terms in <i>Af which vary with N do not vary 
with position, we can establish a correspondance between 
matrix elements of SHO states and optical states for 
those operators which do not involve derivatives. Ex- 
plicitly, if O is an operator (e.g., some function of x), 
then 



(m, n\ O \mn] 



opt 



m,n\0\mn)^^^ . (F7) 



d. Creation and annihilation operators for SHO states 

In quantum mechanics, many computations involving 
simple harmonic oscillator states can be rewritten in a 
more readable form by the use of creation and annihi- 
lation operators |^ 0, 0| . These simplifications occur 
because we can relate any excited state |0) to the ground 
state via the action of these operators: 



/ sho 

at 



|0)sho 



±(l_A^ 
^/2 \A dx 



(F8) 
(F9) 



2. Tabulating useful matrix elements 

When performing perturbation-theory calculations, we 
need the matrix elements of various operators (i.e. x and 
x^) relative to the Hermite-Gauss basis. These elements 
can be expressed as follows: 



= i{N-N)e/2 



(FlOa) 



m, n 



A {cLx + al) 



^/2 



sho 



V2 



'm.m — 1 "l" G 



(m, n X m,n 



opt 



(FlOb) 



= An~N)0/2 



m, n 



A^ + 



sho 



2 ^nn 



e'^ VTO(m+T)(5TO,,yi^2 

+ {2m+ l)5yn..rh 



e '^^y'rhijh + l)5f, 



In the above, we use the subscript opt to denote an 
optical state and sho to denote a SHO state. 



3. Effects of tilt on beam state 

If we tilt mirror 2 about the y axis, the mirror's height 
changes according to 82 = 2kx(f) [cf. Eqs. Ht)l|l and (|30|) ]. 
Using that expression in the two perturbative expansion 
expressions, Eqs. (|32b|l and H32c|l . and taking account of 
the above tools, we can find explicit expressions for the 
first and second order changes in the beam state evalu- 
ated at mirror 1. 



|i,o) (0VZfc)- 



ie/2 



V2 (l-ff2)3/4 



|2,0) {(jyVklf 



(I -.9)^2(1 -.g2) 



(Fll) 



(F12) 



In performing the above calculation, we have explicitly 
substituted in A = A{z±) from Eq. ^Vic\ and e^^ = g + 
i\J\ — [derived from Eq. l)F3b|l ] when appropriate. 



4. Effect of displacement on beam state 

When the mirror is displaced through a distance d 
in the X direction, the mirror height is perturbed from 
h(r) = r2/27^ into 

9 

"XS S 

h'{r) - h{r + sx) = h{r) + — + —. (F13) 

We can in a straightforward manner insert (5/idisp = h' — h 
into the perturbative expansion for the state [Eq. ljE7b() ]: 
note that 6h has terms both first and second order in 
tp; and thereby deduce the appropriate first and second 
order changes in the optical state of the cavity. 

We will not provide a comprehensive calculation here. 
Instead, we note that the first-order changes in height 
6h due to tilt and displacement are identical, involving 
merely a change in factor: 



0Vlfc -> (s/b) (1 - g) 



(F14) 



Therefore, we deduce the first-order effect of a displace- 
ment of mirror 2 is 

= -|l,0)(rf/&) ^ J,L (FI5) 



V2 (l-ff2)3/4 
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APPENDIX G: POSTPROCESSING 
PERFORMED ON FINITE-ELEMENT 
COMPUTATIONS USED IN THERMOELASTIC 
NOISE INTEGRALS 

We used a commercial finite-element code, FEMLAB 
[25|, to solve the fluctuation-dissipation-motivated ax- 
isymmetric elastic model problem presented in Eq. (|f 0|l . 
This code expresses its solution for in terms of a piece- 
wise linear function. Unfortunately, the thermoelastic 
noise integral depends on second derivatives oiya- There- 
fore, even though the FEMLAB program includes many 
useful tools to evaluate physically interesting quantities, 
we had to evaluate the thermoelastic noise integral our- 
selves. We therefore sent the FEMLAB output for y"^ to 
a two-dimensional MATLAB array, then used our own 
MATLAB routines to generate expressions for |V0p and 
thus Ia- 

The two groups performing this computation (Caltech 
and Moscow) each developed relatively independent post- 
processing codes. In this appendix, we describe the Cal- 
tech postprocessing code. 



numbers of points in the r and z directions were 
chosen as ^jNRi/H and ^JNH/Ri, respectively. 

3. Select cutojf region for grid: Because we used a 
grid containing too large a region, we must provide 
a cutoff filter to select only those gridpoints which 
contain physical values. Furthermore, as a practi- 
cal matter, the FEMLAB code gave odd results'^^ 
when we evaluated the response at the edges of the 
computational domain. We therefore chose a cut- 
off filter which eliminated all exterior and "near-to- 
the-boundary" points. 

This filter was applied repeatedly (i.e. during each 
derivative process). Because the filter was such a 
universal ingredient to each subsequent action, we 
will not mention each occurrence on which it is used 
in the following postprocessing. 

4. Compute expansion: Using the relationship be- 
tween 6 and y°', namely 

e^d,V + drU + U/r (Gl) 



1. Caltech postprocessing code 

The Caltech group cobbled together a rough set of 
MATLAB routines to convert the finite-element represen- 
tation of the response to useful values for I a- While 
this code contains more than a few kludges used to get it 
to work, we found the combination of the finite-element 
solution and this postprocessing routine gave answers for 
I A that agreed well with and converged to the values of 
Ia for known, analytic solutions (i.e. the LT cylindrical 
solutions). 



a. Postprocessing technique 



we compute Q. Derivatives are formed as cen- 
tered differences, which are then interpolated (or, 
for the endpoints, extrapolated) back to the grid- 
points. The values of J7/r at r = are found by 
extrapolating the values of U/r for r 7^ 0. 

[To circumvent problems that arose due to dividing 
by small numbers near r = 0, we typically erased 
the values for Q we obtained on four gridlines near 
r — and replaced them by extrapolated values of 
the region immediately outside.] 

5. Compute |V8p: Next, we computed the two 
derivatives drll and d^V and used them to form 



1. Preliminaries: We used the FEMLAB finite- 
element code, with some number N elements, to 
find the elastic response of a frustum of front ra- 
dius back radius R2, and thickness H. We 
increased N until we felt comfortable with the re- 
sulting response curve. The FEMLAB code worked 
in cylindrical coordinates (i.e. r,z, (/)), and assumed 
axisymmetry. 

2. Output to rectangular grid: The code allowed us to 
extract the values of the displacement vector (com- 
ponents U — yr and V — yz) at any point. We 
used this ability to obtain values of U and V on 
a rectangular grid. The code naturally provided 
smoothly extrapolated values of U and V when we 
asked for a point outside the volume simulated. 



|Ve|2 = (drUf + (dzVf . (G2) 

As before, derivatives were evaluated with centered 
differences which were interpolated to gridpoints. 

6. Compute I a: Finally, we used our own two- 
dimensional (Simpson's rule) integrator to evaluate 
the two-dimensional integral 

Ia-- J rdrdz |Vep , (G3) 

which provides an explicit form for Eq. @. 



The total number of points on the rectangular grid 39 i„ ^^^er words, we found NaN (i.e. "not a number") answers 

was chosen to be comparable to N . Specifically, the from our code. 
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b. Testing the result 

Our postprocessing code is far from polished. To in- 
sure the cutoff filter is operating properly and to oth- 
erwise guarantee that the results of the postprocessing 
code seem physically plausible, we usually plotted |V9p 
to verify that the integrand is indeed a well-behaved (i.e. 
smooth-looking) function. 

Because we had a serious limitation on the number 
of points we could practically employ in a reasonable 
amount of desktop computing time (~ few xlO'*), we 
did not perform systematic convergence testing. How- 
ever, what testing we did, corroborated by comparisons 
between our numerical method and exact analytic solu- 
tions, suggests our results are relatively accurate. More 
critically, our computations agreed well with independent 
computations performed by Sergey Strigin and Sergey 
Vyatchanin. 



where H{x) is a step function which is 1 when x > and 
otherwise. 

Because we have complete transverse translation sym- 
metry, we can make our results more tractable by fourier- 
transforming in the transverse dimensions: 

e{K,z) = J e-''^ '^eiR,z)d^R, (H2) 

P{K) = J e-''^ "P{R)(fR . (H3) 

For example, the convolution relating light intensity pro- 
file to the associated elastic response, Eq. (|Hla|l . can be 
re-expressed as 

Q{K,z) = G^{z,K)P{K) (H4) 

&^Uz,k) = _(i±4(i_Me-l-^l . (H5) 



APPENDIX H: THERMOELASTIC NOISE OF 
HALF-INFINITE MIRRORS 

To evaluate the thermoelastic noise associated with a 
given beam shape P{r), we must evaluate the integral I a 
[Eq. (jSJ] given the solution to a model elasticity prob- 
lem [Eq. H10|l ]. As discussed in Sec. lV ASl if the mirror is 
sufficiently large compared to the beam shape P{r), we 
can effectively treat the mirror as half-infinite (i.e. fill- 
ing the whole volume z < 0) in the elasticity problem. 
In this case, the bulk acceleration term in the elasticity 
problem drops out [i.e. Va — > oo in Eq. HlOafl ] and we seek 
only the elastic response of a half-infinite medium to an 
imposed surface stress. This last problem has been dis- 
cussed extensively in the literature — cf., e.g., — 
and there exist simple fourier-based computational tech- 
niques to generate and manage solutions. We apply these 
known solutions from the literature to evaluate the ther- 
moelastic integral Ia- 



1. Elastic solutions for the expansion (Q) 

In the case of half-infinite mirrors, the response to 
the imposed pressure profile P{r) can be found in the 
literature [cf. Eqs. (8.18) and (8.19) of Landau and Lif- 
shitz's book on elasticity where, however, the half- 
infinite volume is chosen above the z = plane rather 
than below; see also Nakagawa et al. , especially their 
Appendix A] . These expressions allow us to explicitly re- 
late the expansion Q to the imposed pressure profile P(r): 



e (f, z) = / G'®^ (r, z; r') P (/) (fP (Hla) 
(1 + cr) (1 - 2cr) zH{-z) 



2. Thermoelastic integral I a 

Inserting the solution discussed above into Eq. © and 
using fourier-transform techniques to simplify the result- 
ing integral, we find 



I A - 



(1 + g) (1 - 2(j) 
2ttE 



d^K \K\ 



P{K) . (H6) 



G(e)( 



27r£;|(f-r;)2 + z2|3/2 



(Hlb) 



APPENDIX I: CONFIGURATIONS WITH 
NONIDENTICAL FRUSTUM MIRRORS AND 
MESA-LIKE BEAMS 

Physically the ETM and ITM need not be identical: 
while the ITM back face size is constrained to be above a 
certain radius (for a given mesa-beam D) by diffraction 
losses, the ETM back face has no size restriction. There- 
fore, if we allow different mirrors, we introduce more pa- 
rameters and therefore more possibilities for finding a 
mirror and beam configuration with very low noise. 

Unfortunately, strictly speaking the mesa beams pre- 
sented in this paper are designed specifically for identical 
cavities. Therefore, in the main text, when we explored 
how thermoelastic noise varied with mirror dimensions 
and shape, we restricted attention to cavities bounded 
by identical mirrors. 

Nonetheless, we expect that generic flat-topped beams 
will always closely resemble the mesa-beam profile, with 
roughly the mesa-beam diffraction losses. Therefore, this 
appendix generalizes the computations of Sec. IVII to in- 
clude distinct ETM and ITM mirrors. 

More specifically, we (i) assume that, for any Di 
and I?2, there is some flat-topped beam whose inten- 
sity very closely mimics the mesa-beam intensity profile 
|umesa(r, £)i)|^ at mirror 1 and |wmesa(?', -D2)|^ at mirror 
2; (ii) select Di, D2, Ri, and R2 so the clipping approxi- 
mation to the diffraction losses is 10 ppm at each mirror 
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face (Sec. Ill D|) : and (iii) use the resulting beam intensity 
profiles to compute the thermoelastic noise integrals Ii 
and h (Sec. ESI- The power spectrum of thermoelastic 
noise for an interferometer with two identical arm cavi- 
ties, each of which has an ITM like mirror 1 and an ETM 
mirror 2 can be expressed relative to the baseline [cf. Eq. 
0] by 



ilTM 



'ETM 



2/bi 



(11) 



1. Thermoelastic noise integral for ETM designs 

The ETM mirror has only two design constraints: (i) 
its mass must be 40 kg and (ii) the diffraction losses off 
its front face must be less than or equal to 10 ppm. Table 
I VII shows the value of the thermoelastic noise integral for 
various values of Rpi, Rp2, and D, chosen to satisfy the 
various constraints (mirror mass 40 kg; diffraction losses 
off the front face of 10 ppm). 

While we have complete freedom to adjust the back 
face size of the ETM, from the table, we see the back 
face size matters little so long as it is less than Rp2 < 9 
cm. By fitting a quadratic to {D/b, I/Ibl) [for only those 
ETM entries with Rp2 ^ 9cm], we estimate the optimal 
ETM dimensions; our results are tabulated in Table IVIll 

In Table En we also combine the optimal ETM with 
the optimal ITM and use Eq. (|lT|l to express the ther- 
moelastic noise of a nonidcntical-mirror configuration rel- 
ative to the baseline. When the ETM can be different 
from the ITM, the optimal thermoelastic noise is signif- 
icantly reduced (because the ETM inner radius can be 
larger) . 



2. Aside: Constructing mesa-like beams 
appropriate to nonidentical mirrors 

In this appendix, we assume that for any Di and D2, 
there is some flat-topped beam whose intensity closely 
mimics the mesa-beam intensity profile |Mniesa(?', -Di)]^ at 
mirror 1 and |umesa('', D2)\^ at mirror 2. In fact, we may 
explicitly construct such a beam through a generalization 
of the mesa-beam proposal (though the authors have not 
numerically explored this possibility). 

The mesa beams constructed in this paper were built 
by averaging minimal gaussian beams which all travelled 
parallel to the optic axis [cf. Eq. |SJ|]. More generally, 
we can construct mesa-like beams by averaging gaussians 
which travel at a position-dependent angle to the optic 
axis; for example, at mirror 1, we could choose the beams 
to have angle distribution 



R2 ^ R^ 



Ri 



(12) 



TABLE VI: The thermoelastic integral / for a frustum end 
test mass (ETM) and a Mexican-Hat Beam, in units of Jbl = 
2.57 X 10"^**s*g"^cm"\ The values of ///bl are estimated 
to be accurate to within one per cent. 



Ri 




-Rpi [cm 


-Rp2 [cm] 


-^[cm] 


D/b 


///bl C 


1 [ppm 


Rpi 


— 8mm 


17.11 


2.00 


28.85 


4.00 


0.175 


10 




— 8mm 


17.11 


8.00 


19.34 


4.00 


0.175 


10 


Rpi 


— 8mm 


17.11 


12.00 


14.87 


4.00 


0.192 


10 


R„^ 


— 8mm 


19.58 


3.00 


21.17 


5.00 


0.134 


10 




— 8mm 


19.58 


9.00 


14.91 


5.00 


0.133 


10 


Rpi 


— 8mm 


19.58 


13.00 


11.83 


5.00 


0.180 


10 


R„^ 


— 8mm 


19.98 


5.00 


18.22 


5.15 


0.133 


10 




— 8mm 


22.08 


3.00 


16.97 


6.00 


0.156 


10 




— 8mm 


22.08 


9.00 


12.45 


6.00 


0.157 


10 


R„^ 


— 8mm 


22.08 


15.00 


9.15 


6.00 


0.310 


10 


Rpi 




13.94 


5.00 


32.93 


3.00 


0.333 


10 


Rpi 




13.94 


7.50 


26.80 


3.00 


0.340 


10 


Rpi 




13.94 


10.00 


21.96 


3.00 


0.345 


10 






16.37 


2.00 


31.34 


4.00 


0.161 


10 


Rpi 




16.37 


4.00 


27.33 


4.00 


0.160 


10 


Rpi 




16.37 


6.00 


23.74 


4.00 


0.160 


10 






16.37 


8.00 


20.63 


4.00 


0.160 


10 


Rpi 




16.37 


10.00 


17.96 


4.00 


0.162 


10 


Rpi 




16.37 


12.00 


15.70 


4.00 


0.173 


10 


Rpi 




16.37 


14.00 


13.78 


4.00 


0.207 


10 


R 




18.85 


1.00 


25.45 


5.00 


0.112 


10 


R„i 

^ •■pi 




18.85 


3.00 


22.69 


5.00 


0.112 


10 






18.85 


5.00 


20.12 


5.00 


0.112 


10 


-Rpi 




18.85 


7.00 


17.81 


5.00 


0.111 


10 






18.85 


9.00 


15.56 


5.00 


0.113 


10 


-Rpi 




18.85 


11.00 


13.97 


5.00 


0.119 


10 


-Rpi 




18.85 


13.00 


12.41 


5.00 


0.144 


10 


jtpi 




19.86 


5.00 


18.42 


5.42 


0.107 


1 n 


-Rpi 




21.36 


3.00 


18.04 


6.00 


0.116 


10 


-Rpi 




21.36 


5.00 


16.24 


6.00 


0.115 


10 


Rpi 




21.36 


7.00 


15.39 


6.00 


0.115 


10 


Rpi 




21.36 


9.00 


14.58 


6.00 


0.116 


10 


Rpi 




21.36 


12.00 


11.15 


6.00 


0.133 


10 


Rpi 




21.36 


15.00 


9.55 


6.00 


0.224 


10 


"D 
(rad 


flfraction losses in each bounce off inner face of the test 
us Ri), in ppm (parts per million). 


mass 



The resulting mesa-type beams will have power dis- 
tributions at the two mirrors which closely resemble the 



power distributions of mesa beams: (i) they will be flat 
over a large central region, of size Di at mirror 1 and size 
D2 ~ Di{R2/ Ri) at mirror 2; and (ii) they will fall of 
rapidly outside this region, with a falloff rate nearly the 
same as that of the minimal gaussian. 
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TABLE VII: Optimized test-mass and light beam configu- 
rations, their thermoelastic noise compared to the baseline. 
[A subset of this table appears as Table I in MBI and as 
Table EVl in this paper.] 
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